Sweetgum Genetic Study

General info

I created this website to track the progress with the sweetgum genetic study.

There are 7 parts:

1. Libraries, dataset loading, and data wrangling.
2. SAS code.
3. Survival rate plots.
4. Groundline diameter (GLD) plots and stats.
5. Tree height plots and stats.
6. DBH plots and stats.
7. Live crown height plots and stats.

Data Wrangling

library(tidyverse)

## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr     1.1.4     ✔ readr     2.1.4
## ✔ forcats   1.0.0     ✔ stringr   1.5.1
## ✔ ggplot2   3.4.4     ✔ tibble    3.2.1
## ✔ lubridate 1.9.3     ✔ tidyr     1.3.0
## ✔ purrr     1.0.2     
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors

library(ggpubr)
library(plotly)

## 
## Attaching package: 'plotly'
## 
## The following object is masked from 'package:ggplot2':
## 
##     last_plot
## 
## The following object is masked from 'package:stats':
## 
##     filter
## 
## The following object is masked from 'package:graphics':
## 
##     layout

library(htmltools)

hane_crown <- read.csv("https://raw.githubusercontent.com/azd169/R-datasets/main/hane_crown.csv",
                       header = T)
hane_gld <- read.csv("https://raw.githubusercontent.com/azd169/R-datasets/main/hane_gld.csv",
                       header = T)
hane_dbh <- read.csv("https://raw.githubusercontent.com/azd169/R-datasets/main/hane_dbh.csv",
                       header = T)
hane_ht <- read.csv("https://raw.githubusercontent.com/azd169/R-datasets/main/hane_ht.csv",
                       header = T)
hane_survival_year <- read.csv("https://raw.githubusercontent.com/azd169/R-datasets/main/hane_survival_year.csv",
                       header = T)
hane_survival_year_family <- read.csv("https://raw.githubusercontent.com/azd169/R-datasets/main/hane_survival_year_family.csv",
                       header = T)
hane_survival_year_location <- read.csv("https://raw.githubusercontent.com/azd169/R-datasets/main/hane_survival_year_location.csv",
                       header = T)
hane_survival_year_location_family <- read.csv("https://raw.githubusercontent.com/azd169/R-datasets/main/hane_survival_year_location_family.csv",
                       header = T)

hane_crown <- hane_crown %>%
  mutate(Year = fct_relevel(Year, c("Fifth", "Sixth", "Seventh")))

hane_gld <- hane_gld %>%
  mutate(Year = fct_relevel(Year, c("Initial", "First", "Second")))

hane_dbh <- hane_dbh %>%
  mutate(Year = fct_relevel(Year, c("Third", "Fourth", "Fifth", "Sixth", "Seventh")))

hane_ht <- hane_ht %>%
  mutate(Year = fct_relevel(Year, c("Initial", "First", "Second", "Third", "Fourth", "Fifth", "Sixth", "Seventh")))

hane_survival_year <- hane_survival_year %>%
  mutate(Year = fct_relevel(Year, c("Initial", "First", "Second")))

hane_survival_year_family <- hane_survival_year_family %>%
  mutate(Year = fct_relevel(Year, c("Initial", "First", "Second")))

hane_survival_year_location <- hane_survival_year_location %>%
  mutate(Year = fct_relevel(Year, c("Initial", "First", "Second")))

hane_survival_year_location_family <- hane_survival_year_location_family %>%
  dplyr::mutate(ID = paste(Location, Year, sep = ", ")) %>%
  mutate(ID = fct_relevel(ID, c("Hill Farm, Initial", "Hill Farm, First", "Hill Farm, Second",
                                "LA Tech, Initial", "LA Tech, First", "LA Tech, Second")))

SAS Code

Procedure:

External trees were removed from the analysis and all the Variables were averaged at plot level. A unique ID was created combining Genotype, Location, and Block (GLD in the SAS code)

Datasets used:

1. WORK.DBH for DBH
2. WORK.GLD for GLD
3. WORK.HEIGHT for tree height
4. WORK.CROWN for crown height

Here are the SAS models:

1. GLD

PROC GLIMMIX DATA = GLD H PLOTS = ALL; CLASS Year Location Block Family GLB; MODEL GLD = Year|Location|Family / SOLUTION DDFM = KR; RANDOM intercept / SOLUTION SUBJECT = GLB TYPE = AR(1) CL; LSMEANS Location Family Year * Location Year * Family/ ADJUST = TUKEY LINES; ODS GRAPHICS OFF; RUN;

2. Height

PROC GLIMMIX DATA = HEIGHT H PLOTS = ALL; CLASS Year Location Block Family GLB; MODEL Height = Year|Location|Family / SOLUTION DDFM = KR; RANDOM intercept / SOLUTION SUBJECT = GLB TYPE = AR(1) CL; LSMEANS Family/ ADJUST = TUKEY LINES; SLICE Year * Location / SLICEBY = Year ADJUST = TUKEY LINES DIFF; SLICE Year * Family / SLICEBY = Year ADJUST = TUKEY LINES DIFF; ODS GRAPHICS OFF; RUN;

3. DBH

PROC GLIMMIX DATA = DBH H PLOTS = ALL; CLASS Year Location Block Family GLB; MODEL DBH = Year|Location|Family / SOLUTION DDFM = KR; RANDOM intercept / SOLUTION SUBJECT = GLB TYPE = AR(1) CL; LSMEANS Family/ ADJUST = TUKEY LINES; SLICE Year * Location / SLICEBY = Year ADJUST = TUKEY LINES DIFF; SLICE Year * Family / SLICEBY = Year ADJUST = TUKEY LINES DIFF; ODS GRAPHICS OFF; RUN;

4. Crown Height

PROC GLIMMIX DATA = CROWN H PLOTS = ALL; CLASS Year Location Block Family GLB; MODEL CROWN = Year|Location|Family / SOLUTION DDFM = KR; RANDOM intercept / SOLUTION SUBJECT = GLB TYPE = AR(1) CL; LSMEANS Family/ ADJUST = TUKEY LINES; SLICE Year * Location / SLICEBY = Year ADJUST = TUKEY LINES DIFF; SLICE Year * Family / SLICEBY = Year ADJUST = TUKEY LINES DIFF; ODS GRAPHICS OFF; RUN;

SAS Output GLD

htmltools::includeHTML("https://raw.githubusercontent.com/azd169/SAS_outputs/main/hane-results-GLD.html")

## Warning: `includeHTML()` was provided a `path` that appears to be a complete HTML document.
## ✖ Path: https://raw.githubusercontent.com/azd169/SAS_outputs/main/hane-results-GLD.html
## ℹ Use `tags$iframe()` to include an HTML document. You can either ensure `path` is accessible in your app or document (see e.g. `shiny::addResourcePath()`) and pass the relative path to the `src` argument. Or you can read the contents of `path` and pass the contents to `srcdoc`.

Results: hane.sas

The GLIMMIX Procedure

The Glimmix Procedure

Model Information

Model Information
Data Set	WORK.GLD
Response Variable	GLD
Response Distribution	Gaussian
Link Function	Identity
Variance Function	Default
Variance Matrix Blocked By	GLB
Estimation Technique	Restricted Maximum Likelihood
Degrees of Freedom Method	Kenward-Roger
Fixed Effects SE Adjustment	Kenward-Roger

Class Level Information

Class Level Information
Class	Levels	Values
Year	3	First Initial Second
Location	2	Hill Farm LA Tech
Block	4	1 2 3 4
Family	6	AGH2 AGH25 AGHS1 AGHS2 AGHS3 AGHS4
GLB	48	AGH25_Hill Farm_1 AGH25_Hill Farm_2 AGH25_Hill Farm_3 AGH25_Hill Farm_4 AGH25_LA Tech_1 AGH25_LA Tech_2 AGH25_LA Tech_3 AGH25_LA Tech_4 AGH2_Hill Farm_1 AGH2_Hill Farm_2 AGH2_Hill Farm_3 AGH2_Hill Farm_4 AGH2_LA Tech_1 AGH2_LA Tech_2 AGH2_LA Tech_3 AGH2_LA Tech_4 AGHS1_Hill Farm_1 AGHS1_Hill Farm_2 AGHS1_Hill Farm_3 AGHS1_Hill Farm_4 AGHS1_LA Tech_1 AGHS1_LA Tech_2 AGHS1_LA Tech_3 AGHS1_LA Tech_4 AGHS2_Hill Farm_1 AGHS2_Hill Farm_2 AGHS2_Hill Farm_3 AGHS2_Hill Farm_4 AGHS2_LA Tech_1 AGHS2_LA Tech_2 AGHS2_LA Tech_3 AGHS2_LA Tech_4 AGHS3_Hill Farm_1 AGHS3_Hill Farm_2 AGHS3_Hill Farm_3 AGHS3_Hill Farm_4 AGHS3_LA Tech_1 AGHS3_LA Tech_2 AGHS3_LA Tech_3 AGHS3_LA Tech_4 AGHS4_Hill Farm_1 AGHS4_Hill Farm_2 AGHS4_Hill Farm_3 AGHS4_Hill Farm_4 AGHS4_LA Tech_1 AGHS4_LA Tech_2 AGHS4_LA Tech_3 AGHS4_LA Tech_4

Number of Observations

Number of Observations Read	144
Number of Observations Used	144

Dimensions

Dimensions
G-side Cov. Parameters	2
R-side Cov. Parameters	1
Columns in X	84
Columns in Z per Subject	1
Subjects (Blocks in V)	48
Max Obs per Subject	3

Optimization Information

Optimization Information
Optimization Technique	Dual Quasi-Newton
Parameters in Optimization	2
Lower Boundaries	2
Upper Boundaries	1
Fixed Effects	Profiled
Residual Variance	Profiled
Starting From	Data

Iteration History

Iteration History
Iteration	Restarts	Evaluations	Objective Function	Change	Max Gradient
0	0	4	640.72210261	.	1.46E-13

Convergence Status

Convergence criterion (ABSGCONV=0.00001) satisfied.

Hessian Matrix

Hessian Matrix
Cov Parm	Subject	Row	Col1	Col2
Variance	GLB	1	80.7379
AR(1)	GLB	2

Fit Statistics

Fit Statistics
-2 Res Log Likelihood	640.72
AIC (smaller is better)	646.72
AICC (smaller is better)	646.95
BIC (smaller is better)	652.34
CAIC (smaller is better)	655.34
HQIC (smaller is better)	648.84
Generalized Chi-Square	1275.07
Gener. Chi-Square / DF	11.81

Covariance Parameter Estimates

Covariance Parameter Estimates
Cov Parm	Subject	Estimate	Standard Error
Variance	GLB	2.5015	1.6529
AR(1)	GLB	0	.
Residual		11.8062	1.9677

Solutions for Fixed Effects

Solutions for Fixed Effects
Effect	Year	Location	Family	Estimate	Standard Error	DF	t Value	Pr > |t|
Intercept				48.3072	1.8913	101.8	25.54	<.0001
Year	First			-34.9458	2.4296	72	-14.38	<.0001
Year	Initial			-43.5014	2.4296	72	-17.90	<.0001
Year	Second			0	.	.	.	.
Location		Hill Farm		-8.7071	2.6747	101.8	-3.26	0.0015
Location		LA Tech		0	.	.	.	.
Year*Location	First	Hill Farm		9.0771	3.4360	72	2.64	0.0101
Year*Location	First	LA Tech		0	.	.	.	.
Year*Location	Initial	Hill Farm		9.6024	3.4360	72	2.79	0.0067
Year*Location	Initial	LA Tech		0	.	.	.	.
Year*Location	Second	Hill Farm		0	.	.	.	.
Year*Location	Second	LA Tech		0	.	.	.	.
Family			AGH2	-12.1855	2.6747	101.8	-4.56	<.0001
Family			AGH25	-6.4656	2.6747	101.8	-2.42	0.0174
Family			AGHS1	-0.8208	2.6747	101.8	-0.31	0.7596
Family			AGHS2	2.9486	2.6747	101.8	1.10	0.2729
Family			AGHS3	-1.5495	2.6747	101.8	-0.58	0.5637
Family			AGHS4	0	.	.	.	.
Year*Family	First		AGH2	9.1584	3.4360	72	2.67	0.0095
Year*Family	First		AGH25	3.8052	3.4360	72	1.11	0.2718
Year*Family	First		AGHS1	-0.1953	3.4360	72	-0.06	0.9548
Year*Family	First		AGHS2	-2.9219	3.4360	72	-0.85	0.3979
Year*Family	First		AGHS3	0.4433	3.4360	72	0.13	0.8977
Year*Family	First		AGHS4	0	.	.	.	.
Year*Family	Initial		AGH2	12.1499	3.4360	72	3.54	0.0007
Year*Family	Initial		AGH25	6.5162	3.4360	72	1.90	0.0619
Year*Family	Initial		AGHS1	0.5714	3.4360	72	0.17	0.8684
Year*Family	Initial		AGHS2	-3.2294	3.4360	72	-0.94	0.3504
Year*Family	Initial		AGHS3	1.6059	3.4360	72	0.47	0.6417
Year*Family	Initial		AGHS4	0	.	.	.	.
Year*Family	Second		AGH2	0	.	.	.	.
Year*Family	Second		AGH25	0	.	.	.	.
Year*Family	Second		AGHS1	0	.	.	.	.
Year*Family	Second		AGHS2	0	.	.	.	.
Year*Family	Second		AGHS3	0	.	.	.	.
Year*Family	Second		AGHS4	0	.	.	.	.
Location*Family		Hill Farm	AGH2	6.2013	3.7826	101.8	1.64	0.1042
Location*Family		Hill Farm	AGH25	1.1020	3.7826	101.8	0.29	0.7714
Location*Family		Hill Farm	AGHS1	2.9179	3.7826	101.8	0.77	0.4423
Location*Family		Hill Farm	AGHS2	-3.6520	3.7826	101.8	-0.97	0.3366
Location*Family		Hill Farm	AGHS3	-1.4682	3.7826	101.8	-0.39	0.6987
Location*Family		Hill Farm	AGHS4	0	.	.	.	.
Location*Family		LA Tech	AGH2	0	.	.	.	.
Location*Family		LA Tech	AGH25	0	.	.	.	.
Location*Family		LA Tech	AGHS1	0	.	.	.	.
Location*Family		LA Tech	AGHS2	0	.	.	.	.
Location*Family		LA Tech	AGHS3	0	.	.	.	.
Location*Family		LA Tech	AGHS4	0	.	.	.	.
Year*Location*Family	First	Hill Farm	AGH2	-4.4378	4.8593	72	-0.91	0.3641
Year*Location*Family	First	Hill Farm	AGH25	1.3168	4.8593	72	0.27	0.7872
Year*Location*Family	First	Hill Farm	AGHS1	0.3735	4.8593	72	0.08	0.9389
Year*Location*Family	First	Hill Farm	AGHS2	5.2515	4.8593	72	1.08	0.2834
Year*Location*Family	First	Hill Farm	AGHS3	3.3401	4.8593	72	0.69	0.4941
Year*Location*Family	First	Hill Farm	AGHS4	0	.	.	.	.
Year*Location*Family	First	LA Tech	AGH2	0	.	.	.	.
Year*Location*Family	First	LA Tech	AGH25	0	.	.	.	.
Year*Location*Family	First	LA Tech	AGHS1	0	.	.	.	.
Year*Location*Family	First	LA Tech	AGHS2	0	.	.	.	.
Year*Location*Family	First	LA Tech	AGHS3	0	.	.	.	.
Year*Location*Family	First	LA Tech	AGHS4	0	.	.	.	.
Year*Location*Family	Initial	Hill Farm	AGH2	-6.4171	4.8593	72	-1.32	0.1908
Year*Location*Family	Initial	Hill Farm	AGH25	-1.7337	4.8593	72	-0.36	0.7223
Year*Location*Family	Initial	Hill Farm	AGHS1	-2.6938	4.8593	72	-0.55	0.5810
Year*Location*Family	Initial	Hill Farm	AGHS2	3.1032	4.8593	72	0.64	0.5251
Year*Location*Family	Initial	Hill Farm	AGHS3	1.0562	4.8593	72	0.22	0.8285
Year*Location*Family	Initial	Hill Farm	AGHS4	0	.	.	.	.
Year*Location*Family	Initial	LA Tech	AGH2	0	.	.	.	.
Year*Location*Family	Initial	LA Tech	AGH25	0	.	.	.	.
Year*Location*Family	Initial	LA Tech	AGHS1	0	.	.	.	.
Year*Location*Family	Initial	LA Tech	AGHS2	0	.	.	.	.
Year*Location*Family	Initial	LA Tech	AGHS3	0	.	.	.	.
Year*Location*Family	Initial	LA Tech	AGHS4	0	.	.	.	.
Year*Location*Family	Second	Hill Farm	AGH2	0	.	.	.	.
Year*Location*Family	Second	Hill Farm	AGH25	0	.	.	.	.
Year*Location*Family	Second	Hill Farm	AGHS1	0	.	.	.	.
Year*Location*Family	Second	Hill Farm	AGHS2	0	.	.	.	.
Year*Location*Family	Second	Hill Farm	AGHS3	0	.	.	.	.
Year*Location*Family	Second	Hill Farm	AGHS4	0	.	.	.	.
Year*Location*Family	Second	LA Tech	AGH2	0	.	.	.	.
Year*Location*Family	Second	LA Tech	AGH25	0	.	.	.	.
Year*Location*Family	Second	LA Tech	AGHS1	0	.	.	.	.
Year*Location*Family	Second	LA Tech	AGHS2	0	.	.	.	.
Year*Location*Family	Second	LA Tech	AGHS3	0	.	.	.	.
Year*Location*Family	Second	LA Tech	AGHS4	0	.	.	.	.

Type III Tests of Fixed Effects

Type III Tests of Fixed Effects
Effect	Num DF	Den DF	F Value	Pr > F
Year	2	72	1477.65	<.0001
Location	1	36	5.24	0.0280
Year*Location	2	72	29.73	<.0001
Family	5	36	3.78	0.0075
Year*Family	10	72	3.10	0.0025
Location*Family	5	36	0.56	0.7268
Year*Location*Family	10	72	0.65	0.7634

Solution for Random Effects

Solution for Random Effects
Effect	Subject	Estimate	Std Err Pred	DF	t Value	Pr > |t|	Alpha	Lower	Upper
Intercept	GLB AGH25_Hill Farm_1	1.0602	1.4398	8.699	0.74	0.4809	0.05	-2.2142	4.3346
Intercept	GLB AGH25_Hill Farm_2	-0.05423	1.4398	8.699	-0.04	0.9708	0.05	-3.3287	3.2202
Intercept	GLB AGH25_Hill Farm_3	-1.0796	1.4398	8.699	-0.75	0.4732	0.05	-4.3540	2.1949
Intercept	GLB AGH25_Hill Farm_4	0.07359	1.4398	8.699	0.05	0.9604	0.05	-3.2008	3.3480
Intercept	GLB AGH25_LA Tech_1	-0.5470	1.4398	8.699	-0.38	0.7131	0.05	-3.8215	2.7274
Intercept	GLB AGH25_LA Tech_2	0.09153	1.4398	8.699	0.06	0.9508	0.05	-3.1829	3.3660
Intercept	GLB AGH25_LA Tech_3	-0.6572	1.4398	8.699	-0.46	0.6593	0.05	-3.9316	2.6172
Intercept	GLB AGH25_LA Tech_4	1.1127	1.4398	8.699	0.77	0.4601	0.05	-2.1617	4.3872
Intercept	GLB AGH2_Hill Farm_1	0.3198	1.4398	8.699	0.22	0.8293	0.05	-2.9546	3.5943
Intercept	GLB AGH2_Hill Farm_2	0.3363	1.4398	8.699	0.23	0.8207	0.05	-2.9381	3.6108
Intercept	GLB AGH2_Hill Farm_3	0.1483	1.4398	8.699	0.10	0.9203	0.05	-3.1261	3.4228
Intercept	GLB AGH2_Hill Farm_4	-0.8045	1.4398	8.699	-0.56	0.5904	0.05	-4.0789	2.4699
Intercept	GLB AGH2_LA Tech_1	0.2711	1.4398	8.699	0.19	0.8550	0.05	-3.0033	3.5456
Intercept	GLB AGH2_LA Tech_2	-0.1484	1.4398	8.699	-0.10	0.9203	0.05	-3.4228	3.1261
Intercept	GLB AGH2_LA Tech_3	0.6397	1.4398	8.699	0.44	0.6677	0.05	-2.6348	3.9141
Intercept	GLB AGH2_LA Tech_4	-0.7625	1.4398	8.699	-0.53	0.6097	0.05	-4.0369	2.5120
Intercept	GLB AGHS1_Hill Farm_1	0.7940	1.4398	8.699	0.55	0.5952	0.05	-2.4805	4.0684
Intercept	GLB AGHS1_Hill Farm_2	-0.4143	1.4398	8.699	-0.29	0.7803	0.05	-3.6888	2.8601
Intercept	GLB AGHS1_Hill Farm_3	0.7606	1.4398	8.699	0.53	0.6105	0.05	-2.5138	4.0350
Intercept	GLB AGHS1_Hill Farm_4	-1.1402	1.4398	8.699	-0.79	0.4495	0.05	-4.4146	2.1342
Intercept	GLB AGHS1_LA Tech_1	-0.6798	1.4398	8.699	-0.47	0.6484	0.05	-3.9543	2.5946
Intercept	GLB AGHS1_LA Tech_2	1.3655	1.4398	8.699	0.95	0.3686	0.05	-1.9090	4.6399
Intercept	GLB AGHS1_LA Tech_3	-0.1045	1.4398	8.699	-0.07	0.9438	0.05	-3.3790	3.1699
Intercept	GLB AGHS1_LA Tech_4	-0.5811	1.4398	8.699	-0.40	0.6963	0.05	-3.8555	2.6933
Intercept	GLB AGHS2_Hill Farm_1	-1.2741	1.4398	8.699	-0.88	0.4000	0.05	-4.5485	2.0003
Intercept	GLB AGHS2_Hill Farm_2	-1.9036	1.4398	8.699	-1.32	0.2199	0.05	-5.1780	1.3709
Intercept	GLB AGHS2_Hill Farm_3	2.0137	1.4398	8.699	1.40	0.1966	0.05	-1.2607	5.2881
Intercept	GLB AGHS2_Hill Farm_4	1.1639	1.4398	8.699	0.81	0.4404	0.05	-2.1105	4.4384
Intercept	GLB AGHS2_LA Tech_1	-1.1906	1.4398	8.699	-0.83	0.4304	0.05	-4.4650	2.0838
Intercept	GLB AGHS2_LA Tech_2	-0.4809	1.4398	8.699	-0.33	0.7463	0.05	-3.7553	2.7935
Intercept	GLB AGHS2_LA Tech_3	0.5704	1.4398	8.699	0.40	0.7015	0.05	-2.7040	3.8448
Intercept	GLB AGHS2_LA Tech_4	1.1011	1.4398	8.699	0.76	0.4647	0.05	-2.1733	4.3755
Intercept	GLB AGHS3_Hill Farm_1	1.5679	1.4398	8.699	1.09	0.3054	0.05	-1.7066	4.8423
Intercept	GLB AGHS3_Hill Farm_2	-1.5387	1.4398	8.699	-1.07	0.3140	0.05	-4.8132	1.7357
Intercept	GLB AGHS3_Hill Farm_3	0.3698	1.4398	8.699	0.26	0.8033	0.05	-2.9046	3.6442
Intercept	GLB AGHS3_Hill Farm_4	-0.3989	1.4398	8.699	-0.28	0.7882	0.05	-3.6734	2.8755
Intercept	GLB AGHS3_LA Tech_1	-0.01348	1.4398	8.699	-0.01	0.9927	0.05	-3.2879	3.2609
Intercept	GLB AGHS3_LA Tech_2	-0.1444	1.4398	8.699	-0.10	0.9224	0.05	-3.4188	3.1301
Intercept	GLB AGHS3_LA Tech_3	-0.3383	1.4398	8.699	-0.23	0.8197	0.05	-3.6127	2.9361
Intercept	GLB AGHS3_LA Tech_4	0.4961	1.4398	8.699	0.34	0.7386	0.05	-2.7783	3.7706
Intercept	GLB AGHS4_Hill Farm_1	-0.4645	1.4398	8.699	-0.32	0.7546	0.05	-3.7389	2.8100
Intercept	GLB AGHS4_Hill Farm_2	1.3958	1.4398	8.699	0.97	0.3585	0.05	-1.8786	4.6702
Intercept	GLB AGHS4_Hill Farm_3	-1.1545	1.4398	8.699	-0.80	0.4440	0.05	-4.4289	2.1199
Intercept	GLB AGHS4_Hill Farm_4	0.2232	1.4398	8.699	0.16	0.8804	0.05	-3.0512	3.4976
Intercept	GLB AGHS4_LA Tech_1	-0.5897	1.4398	8.699	-0.41	0.6920	0.05	-3.8641	2.6847
Intercept	GLB AGHS4_LA Tech_2	0.08132	1.4398	8.699	0.06	0.9562	0.05	-3.1931	3.3558
Intercept	GLB AGHS4_LA Tech_3	0.1598	1.4398	8.699	0.11	0.9142	0.05	-3.1147	3.4342
Intercept	GLB AGHS4_LA Tech_4	0.3486	1.4398	8.699	0.24	0.8143	0.05	-2.9258	3.6230

Slices

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year First	1	101.8	4.04	0.0471

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Location	_Location	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year First	Hill Farm	LA Tech	2.1942	1.0919	101.8	2.01	0.0471	0.0482

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Location	Estimate
Year First	Hill Farm	14.2584	A
Year First
Year First	LA Tech	12.0642	B

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Initial	1	101.8	0.33	0.5645

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Location	_Location	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Initial	Hill Farm	LA Tech	0.6313	1.0919	101.8	0.58	0.5645	0.5650

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Location	Estimate
Year Initial	Hill Farm	5.3606	A
Year Initial			A
Year Initial	LA Tech	4.7294	A

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Second	1	101.8	51.77	<.0001

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Location	_Location	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Second	Hill Farm	LA Tech	-7.8569	1.0919	101.8	-7.20	<.0001	<.0001

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Location	Estimate
Year Second	LA Tech	45.2951	A
Year Second
Year Second	Hill Farm	37.4382	B

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year First	5	101.8	0.77	0.5702

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Family	_Family	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year First	AGH2	AGH25	-0.6944	1.8913	101.8	-0.37	0.7143	0.9991
Year First	AGH2	AGHS1	-2.7749	1.8913	101.8	-1.47	0.1454	0.6859
Year First	AGH2	AGHS2	-2.9718	1.8913	101.8	-1.57	0.1192	0.6198
Year First	AGH2	AGHS3	-1.9752	1.8913	101.8	-1.04	0.2988	0.9012
Year First	AGH2	AGHS4	-2.1454	1.8913	101.8	-1.13	0.2593	0.8654
Year First	AGH25	AGHS1	-2.0806	1.8913	101.8	-1.10	0.2739	0.8798
Year First	AGH25	AGHS2	-2.2774	1.8913	101.8	-1.20	0.2313	0.8334
Year First	AGH25	AGHS3	-1.2808	1.8913	101.8	-0.68	0.4998	0.9839
Year First	AGH25	AGHS4	-1.4510	1.8913	101.8	-0.77	0.4447	0.9721
Year First	AGHS1	AGHS2	-0.1969	1.8913	101.8	-0.10	0.9173	1.0000
Year First	AGHS1	AGHS3	0.7998	1.8913	101.8	0.42	0.6733	0.9982
Year First	AGHS1	AGHS4	0.6296	1.8913	101.8	0.33	0.7399	0.9994
Year First	AGHS2	AGHS3	0.9966	1.8913	101.8	0.53	0.5994	0.9949
Year First	AGHS2	AGHS4	0.8264	1.8913	101.8	0.44	0.6631	0.9979
Year First	AGHS3	AGHS4	-0.1702	1.8913	101.8	-0.09	0.9285	1.0000

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Family	Estimate
Year First	AGHS2	14.3728	A
Year First			A
Year First	AGHS1	14.1760	A
Year First			A
Year First	AGHS4	13.5464	A
Year First			A
Year First	AGHS3	13.3762	A
Year First			A
Year First	AGH25	12.0954	A
Year First			A
Year First	AGH2	11.4010	A

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Initial	5	101.8	0.02	0.9998

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Family	_Family	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Initial	AGH2	AGH25	0.1218	1.8913	101.8	0.06	0.9488	1.0000
Year Initial	AGH2	AGHS1	-0.00609	1.8913	101.8	-0.00	0.9974	1.0000
Year Initial	AGH2	AGHS2	0.4117	1.8913	101.8	0.22	0.8281	0.9999
Year Initial	AGH2	AGHS3	0.006111	1.8913	101.8	0.00	0.9974	1.0000
Year Initial	AGH2	AGHS4	-0.1435	1.8913	101.8	-0.08	0.9397	1.0000
Year Initial	AGH25	AGHS1	-0.1279	1.8913	101.8	-0.07	0.9462	1.0000
Year Initial	AGH25	AGHS2	0.2899	1.8913	101.8	0.15	0.8785	1.0000
Year Initial	AGH25	AGHS3	-0.1157	1.8913	101.8	-0.06	0.9513	1.0000
Year Initial	AGH25	AGHS4	-0.2653	1.8913	101.8	-0.14	0.8887	1.0000
Year Initial	AGHS1	AGHS2	0.4178	1.8913	101.8	0.22	0.8256	0.9999
Year Initial	AGHS1	AGHS3	0.01220	1.8913	101.8	0.01	0.9949	1.0000
Year Initial	AGHS1	AGHS4	-0.1374	1.8913	101.8	-0.07	0.9422	1.0000
Year Initial	AGHS2	AGHS3	-0.4056	1.8913	101.8	-0.21	0.8306	0.9999
Year Initial	AGHS2	AGHS4	-0.5552	1.8913	101.8	-0.29	0.7697	0.9997
Year Initial	AGHS3	AGHS4	-0.1496	1.8913	101.8	-0.08	0.9371	1.0000

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Family	Estimate
Year Initial	AGHS4	5.2535	A
Year Initial			A
Year Initial	AGHS1	5.1161	A
Year Initial			A
Year Initial	AGH2	5.1100	A
Year Initial			A
Year Initial	AGHS3	5.1039	A
Year Initial			A
Year Initial	AGH25	4.9882	A
Year Initial			A
Year Initial	AGHS2	4.6983	A

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Second	5	101.8	9.42	<.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Family	_Family	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Second	AGH2	AGH25	-3.1702	1.8913	101.8	-1.68	0.0968	0.5517
Year Second	AGH2	AGHS1	-9.7229	1.8913	101.8	-5.14	<.0001	<.0001
Year Second	AGH2	AGHS2	-10.2074	1.8913	101.8	-5.40	<.0001	<.0001
Year Second	AGH2	AGHS3	-6.8012	1.8913	101.8	-3.60	0.0005	0.0075
Year Second	AGH2	AGHS4	-9.0848	1.8913	101.8	-4.80	<.0001	0.0001
Year Second	AGH25	AGHS1	-6.5527	1.8913	101.8	-3.46	0.0008	0.0112
Year Second	AGH25	AGHS2	-7.0372	1.8913	101.8	-3.72	0.0003	0.0051
Year Second	AGH25	AGHS3	-3.6310	1.8913	101.8	-1.92	0.0577	0.3987
Year Second	AGH25	AGHS4	-5.9146	1.8913	101.8	-3.13	0.0023	0.0294
Year Second	AGHS1	AGHS2	-0.4845	1.8913	101.8	-0.26	0.7983	0.9998
Year Second	AGHS1	AGHS3	2.9217	1.8913	101.8	1.54	0.1255	0.6368
Year Second	AGHS1	AGHS4	0.6381	1.8913	101.8	0.34	0.7365	0.9994
Year Second	AGHS2	AGHS3	3.4062	1.8913	101.8	1.80	0.0747	0.4716
Year Second	AGHS2	AGHS4	1.1226	1.8913	101.8	0.59	0.5541	0.9912
Year Second	AGHS3	AGHS4	-2.2836	1.8913	101.8	-1.21	0.2301	0.8319

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Family	Estimate
Year Second	AGHS2	45.0763		A
Year Second				A
Year Second	AGHS1	44.5918		A
Year Second				A
Year Second	AGHS4	43.9537		A
Year Second				A
Year Second	AGHS3	41.6701	B	A
Year Second			B
Year Second	AGH25	38.0391	B	C
Year Second				C
Year Second	AGH2	34.8689		C

LS-Means Results

Family LS-Means

Family Least Squares Means
Family	Estimate	Standard Error	DF	t Value	Pr > |t|
AGH2	17.1266	0.8970	36	19.09	<.0001
AGH25	18.3742	0.8970	36	20.48	<.0001
AGHS1	21.2946	0.8970	36	23.74	<.0001
AGHS2	21.3825	0.8970	36	23.84	<.0001
AGHS3	20.0501	0.8970	36	22.35	<.0001
AGHS4	20.9179	0.8970	36	23.32	<.0001

Family Diffs

Differences of Family Least Squares Means Adjustment for Multiple Comparisons: Tukey
Family	_Family	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
AGH2	AGH25	-1.2476	1.2686	36	-0.98	0.3319	0.9202
AGH2	AGHS1	-4.1680	1.2686	36	-3.29	0.0023	0.0255
AGH2	AGHS2	-4.2558	1.2686	36	-3.35	0.0019	0.0214
AGH2	AGHS3	-2.9234	1.2686	36	-2.30	0.0271	0.2187
AGH2	AGHS4	-3.7912	1.2686	36	-2.99	0.0050	0.0524
AGH25	AGHS1	-2.9204	1.2686	36	-2.30	0.0272	0.2196
AGH25	AGHS2	-3.0082	1.2686	36	-2.37	0.0232	0.1934
AGH25	AGHS3	-1.6758	1.2686	36	-1.32	0.1948	0.7716
AGH25	AGHS4	-2.5436	1.2686	36	-2.01	0.0525	0.3594
AGHS1	AGHS2	-0.08785	1.2686	36	-0.07	0.9452	1.0000
AGHS1	AGHS3	1.2445	1.2686	36	0.98	0.3331	0.9210
AGHS1	AGHS4	0.3768	1.2686	36	0.30	0.7682	0.9997
AGHS2	AGHS3	1.3324	1.2686	36	1.05	0.3006	0.8974
AGHS2	AGHS4	0.4646	1.2686	36	0.37	0.7163	0.9991
AGHS3	AGHS4	-0.8678	1.2686	36	-0.68	0.4983	0.9826

Family Lines

Tukey Grouping for Family Least Squares Means (Alpha=0.05)
LS-means with the same letter are not significantly different.
Family	Estimate
AGHS2	21.3825		A
			A
AGHS1	21.2946		A
			A
AGHS4	20.9179	B	A
		B	A
AGHS3	20.0501	B	A
		B	A
AGH25	18.3742	B	A
		B
AGH2	17.1266	B

SAS Output DBH

htmltools::includeHTML("https://raw.githubusercontent.com/azd169/SAS_outputs/main/hane-results-DBH.html")

## Warning: `includeHTML()` was provided a `path` that appears to be a complete HTML document.
## ✖ Path: https://raw.githubusercontent.com/azd169/SAS_outputs/main/hane-results-DBH.html
## ℹ Use `tags$iframe()` to include an HTML document. You can either ensure `path` is accessible in your app or document (see e.g. `shiny::addResourcePath()`) and pass the relative path to the `src` argument. Or you can read the contents of `path` and pass the contents to `srcdoc`.

Results: hane.sas

The GLIMMIX Procedure

The Glimmix Procedure

Model Information

Model Information
Data Set	WORK.DBH
Response Variable	DBH
Response Distribution	Gaussian
Link Function	Identity
Variance Function	Default
Variance Matrix Blocked By	GLB
Estimation Technique	Restricted Maximum Likelihood
Degrees of Freedom Method	Kenward-Roger
Fixed Effects SE Adjustment	Kenward-Roger

Class Level Information

Class Level Information
Class	Levels	Values
Year	5	Fifth Fourth Seventh Sixth Third
Location	2	Hill Farm LA Tech
Block	4	1 2 3 4
Family	6	AGH2 AGH25 AGHS1 AGHS2 AGHS3 AGHS4
GLB	48	AGH25_Hill Farm_1 AGH25_Hill Farm_2 AGH25_Hill Farm_3 AGH25_Hill Farm_4 AGH25_LA Tech_1 AGH25_LA Tech_2 AGH25_LA Tech_3 AGH25_LA Tech_4 AGH2_Hill Farm_1 AGH2_Hill Farm_2 AGH2_Hill Farm_3 AGH2_Hill Farm_4 AGH2_LA Tech_1 AGH2_LA Tech_2 AGH2_LA Tech_3 AGH2_LA Tech_4 AGHS1_Hill Farm_1 AGHS1_Hill Farm_2 AGHS1_Hill Farm_3 AGHS1_Hill Farm_4 AGHS1_LA Tech_1 AGHS1_LA Tech_2 AGHS1_LA Tech_3 AGHS1_LA Tech_4 AGHS2_Hill Farm_1 AGHS2_Hill Farm_2 AGHS2_Hill Farm_3 AGHS2_Hill Farm_4 AGHS2_LA Tech_1 AGHS2_LA Tech_2 AGHS2_LA Tech_3 AGHS2_LA Tech_4 AGHS3_Hill Farm_1 AGHS3_Hill Farm_2 AGHS3_Hill Farm_3 AGHS3_Hill Farm_4 AGHS3_LA Tech_1 AGHS3_LA Tech_2 AGHS3_LA Tech_3 AGHS3_LA Tech_4 AGHS4_Hill Farm_1 AGHS4_Hill Farm_2 AGHS4_Hill Farm_3 AGHS4_Hill Farm_4 AGHS4_LA Tech_1 AGHS4_LA Tech_2 AGHS4_LA Tech_3 AGHS4_LA Tech_4

Number of Observations

Number of Observations Read	240
Number of Observations Used	240

Dimensions

Dimensions
G-side Cov. Parameters	2
R-side Cov. Parameters	1
Columns in X	126
Columns in Z per Subject	1
Subjects (Blocks in V)	48
Max Obs per Subject	5

Optimization Information

Optimization Information
Optimization Technique	Dual Quasi-Newton
Parameters in Optimization	2
Lower Boundaries	2
Upper Boundaries	1
Fixed Effects	Profiled
Residual Variance	Profiled
Starting From	Data

Iteration History

Iteration History
Iteration	Restarts	Evaluations	Objective Function	Change	Max Gradient
0	0	4	337.79560627	.	5.28E-13

Convergence Status

Convergence criterion (ABSGCONV=0.00001) satisfied.

Hessian Matrix

Hessian Matrix
Cov Parm	Subject	Row	Col1	Col2
Variance	GLB	1	7.5640
AR(1)	GLB	2

Fit Statistics

Fit Statistics
-2 Res Log Likelihood	337.80
AIC (smaller is better)	343.80
AICC (smaller is better)	343.93
BIC (smaller is better)	349.41
CAIC (smaller is better)	352.41
HQIC (smaller is better)	345.92
Generalized Chi-Square	27.50
Gener. Chi-Square / DF	0.15

Covariance Parameter Estimates

Covariance Parameter Estimates
Cov Parm	Subject	Estimate	Standard Error
Variance	GLB	0.2675	0.07035
AR(1)	GLB	0	.
Residual		0.1528	0.01800

Solutions for Fixed Effects

Solutions for Fixed Effects
Effect	Year	Location	Family	Estimate	Standard Error	DF	t Value	Pr > |t|
Intercept				3.8518	0.3241	68.68	11.88	<.0001
Year	Fifth			6.8444	0.2764	144	24.77	<.0001
Year	Fourth			3.7379	0.2764	144	13.53	<.0001
Year	Seventh			9.9398	0.2764	144	35.97	<.0001
Year	Sixth			8.7141	0.2764	144	31.53	<.0001
Year	Third			0	.	.	.	.
Location		Hill Farm		-0.9135	0.4584	68.68	-1.99	0.0502
Location		LA Tech		0	.	.	.	.
Year*Location	Fifth	Hill Farm		-1.2948	0.3908	144	-3.31	0.0012
Year*Location	Fifth	LA Tech		0	.	.	.	.
Year*Location	Fourth	Hill Farm		-0.5688	0.3908	144	-1.46	0.1477
Year*Location	Fourth	LA Tech		0	.	.	.	.
Year*Location	Seventh	Hill Farm		-0.6350	0.3908	144	-1.62	0.1064
Year*Location	Seventh	LA Tech		0	.	.	.	.
Year*Location	Sixth	Hill Farm		-0.4200	0.3908	144	-1.07	0.2844
Year*Location	Sixth	LA Tech		0	.	.	.	.
Year*Location	Third	Hill Farm		0	.	.	.	.
Year*Location	Third	LA Tech		0	.	.	.	.
Family			AGH2	-2.1578	0.4584	68.68	-4.71	<.0001
Family			AGH25	-1.9694	0.4584	68.68	-4.30	<.0001
Family			AGHS1	-0.3172	0.4584	68.68	-0.69	0.4913
Family			AGHS2	0.1063	0.4584	68.68	0.23	0.8173
Family			AGHS3	-1.2254	0.4584	68.68	-2.67	0.0094
Family			AGHS4	0	.	.	.	.
Year*Family	Fifth		AGH2	-0.5357	0.3908	144	-1.37	0.1726
Year*Family	Fifth		AGH25	-0.2856	0.3908	144	-0.73	0.4661
Year*Family	Fifth		AGHS1	-1.4184	0.3908	144	-3.63	0.0004
Year*Family	Fifth		AGHS2	-1.6092	0.3908	144	-4.12	<.0001
Year*Family	Fifth		AGHS3	-0.1346	0.3908	144	-0.34	0.7311
Year*Family	Fifth		AGHS4	0	.	.	.	.
Year*Family	Fourth		AGH2	-0.3221	0.3908	144	-0.82	0.4112
Year*Family	Fourth		AGH25	-0.01489	0.3908	144	-0.04	0.9697
Year*Family	Fourth		AGHS1	-0.6544	0.3908	144	-1.67	0.0962
Year*Family	Fourth		AGHS2	-0.7991	0.3908	144	-2.04	0.0427
Year*Family	Fourth		AGHS3	0.06217	0.3908	144	0.16	0.8738
Year*Family	Fourth		AGHS4	0	.	.	.	.
Year*Family	Seventh		AGH2	-0.04498	0.3908	144	-0.12	0.9085
Year*Family	Seventh		AGH25	-0.2542	0.3908	144	-0.65	0.5165
Year*Family	Seventh		AGHS1	-1.5222	0.3908	144	-3.89	0.0002
Year*Family	Seventh		AGHS2	-1.6619	0.3908	144	-4.25	<.0001
Year*Family	Seventh		AGHS3	-0.1730	0.3908	144	-0.44	0.6586
Year*Family	Seventh		AGHS4	0	.	.	.	.
Year*Family	Sixth		AGH2	-0.4924	0.3908	144	-1.26	0.2097
Year*Family	Sixth		AGH25	-0.2868	0.3908	144	-0.73	0.4643
Year*Family	Sixth		AGHS1	-1.2209	0.3908	144	-3.12	0.0022
Year*Family	Sixth		AGHS2	-1.6868	0.3908	144	-4.32	<.0001
Year*Family	Sixth		AGHS3	0.09563	0.3908	144	0.24	0.8070
Year*Family	Sixth		AGHS4	0	.	.	.	.
Year*Family	Third		AGH2	0	.	.	.	.
Year*Family	Third		AGH25	0	.	.	.	.
Year*Family	Third		AGHS1	0	.	.	.	.
Year*Family	Third		AGHS2	0	.	.	.	.
Year*Family	Third		AGHS3	0	.	.	.	.
Year*Family	Third		AGHS4	0	.	.	.	.
Location*Family		Hill Farm	AGH2	0.7712	0.6483	68.68	1.19	0.2383
Location*Family		Hill Farm	AGH25	0.5841	0.6483	68.68	0.90	0.3708
Location*Family		Hill Farm	AGHS1	0.2693	0.6483	68.68	0.42	0.6792
Location*Family		Hill Farm	AGHS2	-0.4456	0.6483	68.68	-0.69	0.4941
Location*Family		Hill Farm	AGHS3	0.2826	0.6483	68.68	0.44	0.6642
Location*Family		Hill Farm	AGHS4	0	.	.	.	.
Location*Family		LA Tech	AGH2	0	.	.	.	.
Location*Family		LA Tech	AGH25	0	.	.	.	.
Location*Family		LA Tech	AGHS1	0	.	.	.	.
Location*Family		LA Tech	AGHS2	0	.	.	.	.
Location*Family		LA Tech	AGHS3	0	.	.	.	.
Location*Family		LA Tech	AGHS4	0	.	.	.	.
Year*Location*Family	Fifth	Hill Farm	AGH2	-0.4343	0.5527	144	-0.79	0.4333
Year*Location*Family	Fifth	Hill Farm	AGH25	-0.2834	0.5527	144	-0.51	0.6089
Year*Location*Family	Fifth	Hill Farm	AGHS1	0.06586	0.5527	144	0.12	0.9053
Year*Location*Family	Fifth	Hill Farm	AGHS2	0.7027	0.5527	144	1.27	0.2057
Year*Location*Family	Fifth	Hill Farm	AGHS3	-0.1706	0.5527	144	-0.31	0.7580
Year*Location*Family	Fifth	Hill Farm	AGHS4	0	.	.	.	.
Year*Location*Family	Fifth	LA Tech	AGH2	0	.	.	.	.
Year*Location*Family	Fifth	LA Tech	AGH25	0	.	.	.	.
Year*Location*Family	Fifth	LA Tech	AGHS1	0	.	.	.	.
Year*Location*Family	Fifth	LA Tech	AGHS2	0	.	.	.	.
Year*Location*Family	Fifth	LA Tech	AGHS3	0	.	.	.	.
Year*Location*Family	Fifth	LA Tech	AGHS4	0	.	.	.	.
Year*Location*Family	Fourth	Hill Farm	AGH2	-0.9294	0.5527	144	-1.68	0.0948
Year*Location*Family	Fourth	Hill Farm	AGH25	-1.0055	0.5527	144	-1.82	0.0710
Year*Location*Family	Fourth	Hill Farm	AGHS1	-0.2051	0.5527	144	-0.37	0.7112
Year*Location*Family	Fourth	Hill Farm	AGHS2	-0.04451	0.5527	144	-0.08	0.9359
Year*Location*Family	Fourth	Hill Farm	AGHS3	-0.7046	0.5527	144	-1.27	0.2044
Year*Location*Family	Fourth	Hill Farm	AGHS4	0	.	.	.	.
Year*Location*Family	Fourth	LA Tech	AGH2	0	.	.	.	.
Year*Location*Family	Fourth	LA Tech	AGH25	0	.	.	.	.
Year*Location*Family	Fourth	LA Tech	AGHS1	0	.	.	.	.
Year*Location*Family	Fourth	LA Tech	AGHS2	0	.	.	.	.
Year*Location*Family	Fourth	LA Tech	AGHS3	0	.	.	.	.
Year*Location*Family	Fourth	LA Tech	AGHS4	0	.	.	.	.
Year*Location*Family	Seventh	Hill Farm	AGH2	-0.8942	0.5527	144	-1.62	0.1079
Year*Location*Family	Seventh	Hill Farm	AGH25	-0.6181	0.5527	144	-1.12	0.2653
Year*Location*Family	Seventh	Hill Farm	AGHS1	-0.3015	0.5527	144	-0.55	0.5862
Year*Location*Family	Seventh	Hill Farm	AGHS2	0.4281	0.5527	144	0.77	0.4398
Year*Location*Family	Seventh	Hill Farm	AGHS3	-0.8381	0.5527	144	-1.52	0.1316
Year*Location*Family	Seventh	Hill Farm	AGHS4	0	.	.	.	.
Year*Location*Family	Seventh	LA Tech	AGH2	0	.	.	.	.
Year*Location*Family	Seventh	LA Tech	AGH25	0	.	.	.	.
Year*Location*Family	Seventh	LA Tech	AGHS1	0	.	.	.	.
Year*Location*Family	Seventh	LA Tech	AGHS2	0	.	.	.	.
Year*Location*Family	Seventh	LA Tech	AGHS3	0	.	.	.	.
Year*Location*Family	Seventh	LA Tech	AGHS4	0	.	.	.	.
Year*Location*Family	Sixth	Hill Farm	AGH2	-0.6702	0.5527	144	-1.21	0.2273
Year*Location*Family	Sixth	Hill Farm	AGH25	-0.7039	0.5527	144	-1.27	0.2049
Year*Location*Family	Sixth	Hill Farm	AGHS1	-0.5868	0.5527	144	-1.06	0.2902
Year*Location*Family	Sixth	Hill Farm	AGHS2	0.3045	0.5527	144	0.55	0.5825
Year*Location*Family	Sixth	Hill Farm	AGHS3	-0.5005	0.5527	144	-0.91	0.3667
Year*Location*Family	Sixth	Hill Farm	AGHS4	0	.	.	.	.
Year*Location*Family	Sixth	LA Tech	AGH2	0	.	.	.	.
Year*Location*Family	Sixth	LA Tech	AGH25	0	.	.	.	.
Year*Location*Family	Sixth	LA Tech	AGHS1	0	.	.	.	.
Year*Location*Family	Sixth	LA Tech	AGHS2	0	.	.	.	.
Year*Location*Family	Sixth	LA Tech	AGHS3	0	.	.	.	.
Year*Location*Family	Sixth	LA Tech	AGHS4	0	.	.	.	.
Year*Location*Family	Third	Hill Farm	AGH2	0	.	.	.	.
Year*Location*Family	Third	Hill Farm	AGH25	0	.	.	.	.
Year*Location*Family	Third	Hill Farm	AGHS1	0	.	.	.	.
Year*Location*Family	Third	Hill Farm	AGHS2	0	.	.	.	.
Year*Location*Family	Third	Hill Farm	AGHS3	0	.	.	.	.
Year*Location*Family	Third	Hill Farm	AGHS4	0	.	.	.	.
Year*Location*Family	Third	LA Tech	AGH2	0	.	.	.	.
Year*Location*Family	Third	LA Tech	AGH25	0	.	.	.	.
Year*Location*Family	Third	LA Tech	AGHS1	0	.	.	.	.
Year*Location*Family	Third	LA Tech	AGHS2	0	.	.	.	.
Year*Location*Family	Third	LA Tech	AGHS3	0	.	.	.	.
Year*Location*Family	Third	LA Tech	AGHS4	0	.	.	.	.

Type III Tests of Fixed Effects

Type III Tests of Fixed Effects
Effect	Num DF	Den DF	F Value	Pr > F
Year	4	144	4058.38	<.0001
Location	1	36	90.58	<.0001
Year*Location	4	144	19.76	<.0001
Family	5	36	18.48	<.0001
Year*Family	20	144	4.79	<.0001
Location*Family	5	36	0.13	0.9849
Year*Location*Family	20	144	0.77	0.7482

Solution for Random Effects

Solution for Random Effects
Effect	Subject	Estimate	Std Err Pred	DF	t Value	Pr > |t|	Alpha	Lower	Upper
Intercept	GLB AGH25_Hill Farm_1	0.2413	0.2963	46.37	0.81	0.4195	0.05	-0.3549	0.8375
Intercept	GLB AGH25_Hill Farm_2	-0.2756	0.2963	46.37	-0.93	0.3570	0.05	-0.8718	0.3206
Intercept	GLB AGH25_Hill Farm_3	-0.01982	0.2963	46.37	-0.07	0.9470	0.05	-0.6160	0.5764
Intercept	GLB AGH25_Hill Farm_4	0.05412	0.2963	46.37	0.18	0.8558	0.05	-0.5421	0.6503
Intercept	GLB AGH25_LA Tech_1	-0.4092	0.2963	46.37	-1.38	0.1738	0.05	-1.0054	0.1870
Intercept	GLB AGH25_LA Tech_2	-0.3500	0.2963	46.37	-1.18	0.2434	0.05	-0.9463	0.2462
Intercept	GLB AGH25_LA Tech_3	0.3857	0.2963	46.37	1.30	0.1994	0.05	-0.2105	0.9819
Intercept	GLB AGH25_LA Tech_4	0.3735	0.2963	46.37	1.26	0.2137	0.05	-0.2227	0.9698
Intercept	GLB AGH2_Hill Farm_1	-0.2716	0.2963	46.37	-0.92	0.3640	0.05	-0.8678	0.3246
Intercept	GLB AGH2_Hill Farm_2	0.1398	0.2963	46.37	0.47	0.6392	0.05	-0.4564	0.7360
Intercept	GLB AGH2_Hill Farm_3	0.3201	0.2963	46.37	1.08	0.2855	0.05	-0.2761	0.9164
Intercept	GLB AGH2_Hill Farm_4	-0.1884	0.2963	46.37	-0.64	0.5280	0.05	-0.7846	0.4079
Intercept	GLB AGH2_LA Tech_1	-0.05234	0.2963	46.37	-0.18	0.8606	0.05	-0.6485	0.5439
Intercept	GLB AGH2_LA Tech_2	-0.6627	0.2963	46.37	-2.24	0.0301	0.05	-1.2590	-0.06653
Intercept	GLB AGH2_LA Tech_3	0.8727	0.2963	46.37	2.95	0.0050	0.05	0.2764	1.4689
Intercept	GLB AGH2_LA Tech_4	-0.1576	0.2963	46.37	-0.53	0.5973	0.05	-0.7538	0.4386
Intercept	GLB AGHS1_Hill Farm_1	0.1936	0.2963	46.37	0.65	0.5166	0.05	-0.4026	0.7898
Intercept	GLB AGHS1_Hill Farm_2	-0.1444	0.2963	46.37	-0.49	0.6283	0.05	-0.7406	0.4518
Intercept	GLB AGHS1_Hill Farm_3	0.4937	0.2963	46.37	1.67	0.1024	0.05	-0.1025	1.0899
Intercept	GLB AGHS1_Hill Farm_4	-0.5429	0.2963	46.37	-1.83	0.0733	0.05	-1.1391	0.05331
Intercept	GLB AGHS1_LA Tech_1	-0.4634	0.2963	46.37	-1.56	0.1246	0.05	-1.0596	0.1328
Intercept	GLB AGHS1_LA Tech_2	0.4700	0.2963	46.37	1.59	0.1195	0.05	-0.1262	1.0662
Intercept	GLB AGHS1_LA Tech_3	0.4483	0.2963	46.37	1.51	0.1370	0.05	-0.1479	1.0445
Intercept	GLB AGHS1_LA Tech_4	-0.4549	0.2963	46.37	-1.54	0.1315	0.05	-1.0511	0.1413
Intercept	GLB AGHS2_Hill Farm_1	-0.00318	0.2963	46.37	-0.01	0.9915	0.05	-0.5994	0.5930
Intercept	GLB AGHS2_Hill Farm_2	-0.8445	0.2963	46.37	-2.85	0.0065	0.05	-1.4407	-0.2483
Intercept	GLB AGHS2_Hill Farm_3	0.3695	0.2963	46.37	1.25	0.2185	0.05	-0.2267	0.9657
Intercept	GLB AGHS2_Hill Farm_4	0.4782	0.2963	46.37	1.61	0.1133	0.05	-0.1180	1.0744
Intercept	GLB AGHS2_LA Tech_1	-0.5213	0.2963	46.37	-1.76	0.0851	0.05	-1.1175	0.07492
Intercept	GLB AGHS2_LA Tech_2	0.05283	0.2963	46.37	0.18	0.8592	0.05	-0.5434	0.6490
Intercept	GLB AGHS2_LA Tech_3	0.5827	0.2963	46.37	1.97	0.0552	0.05	-0.01349	1.1789
Intercept	GLB AGHS2_LA Tech_4	-0.1143	0.2963	46.37	-0.39	0.7015	0.05	-0.7105	0.4820
Intercept	GLB AGHS3_Hill Farm_1	0.8690	0.2963	46.37	2.93	0.0052	0.05	0.2727	1.4652
Intercept	GLB AGHS3_Hill Farm_2	-1.1000	0.2963	46.37	-3.71	0.0005	0.05	-1.6962	-0.5038
Intercept	GLB AGHS3_Hill Farm_3	0.5274	0.2963	46.37	1.78	0.0816	0.05	-0.06881	1.1236
Intercept	GLB AGHS3_Hill Farm_4	-0.2963	0.2963	46.37	-1.00	0.3224	0.05	-0.8926	0.2999
Intercept	GLB AGHS3_LA Tech_1	-0.2220	0.2963	46.37	-0.75	0.4574	0.05	-0.8182	0.3742
Intercept	GLB AGHS3_LA Tech_2	-0.06415	0.2963	46.37	-0.22	0.8295	0.05	-0.6604	0.5321
Intercept	GLB AGHS3_LA Tech_3	0.4623	0.2963	46.37	1.56	0.1254	0.05	-0.1339	1.0585
Intercept	GLB AGHS3_LA Tech_4	-0.1761	0.2963	46.37	-0.59	0.5550	0.05	-0.7724	0.4201
Intercept	GLB AGHS4_Hill Farm_1	0.002375	0.2963	46.37	0.01	0.9936	0.05	-0.5938	0.5986
Intercept	GLB AGHS4_Hill Farm_2	0.5383	0.2963	46.37	1.82	0.0757	0.05	-0.05786	1.1346
Intercept	GLB AGHS4_Hill Farm_3	-0.2894	0.2963	46.37	-0.98	0.3337	0.05	-0.8856	0.3068
Intercept	GLB AGHS4_Hill Farm_4	-0.2513	0.2963	46.37	-0.85	0.4006	0.05	-0.8475	0.3449
Intercept	GLB AGHS4_LA Tech_1	-0.2461	0.2963	46.37	-0.83	0.4104	0.05	-0.8423	0.3501
Intercept	GLB AGHS4_LA Tech_2	-0.2476	0.2963	46.37	-0.84	0.4076	0.05	-0.8438	0.3486
Intercept	GLB AGHS4_LA Tech_3	0.1678	0.2963	46.37	0.57	0.5738	0.05	-0.4284	0.7641
Intercept	GLB AGHS4_LA Tech_4	0.3258	0.2963	46.37	1.10	0.2771	0.05	-0.2704	0.9221

Slices

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Fifth	1	68.68	112.48	<.0001

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Location	_Location	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Fifth	Hill Farm	LA Tech	-1.9848	0.1871	68.68	-10.61	<.0001	<.0001

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Location	Estimate
Year Fifth	LA Tech	9.1051	A
Year Fifth
Year Fifth	Hill Farm	7.1203	B

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Fourth	1	68.68	84.50	<.0001

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Location	_Location	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Fourth	Hill Farm	LA Tech	-1.7203	0.1871	68.68	-9.19	<.0001	<.0001

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Location	Estimate
Year Fourth	LA Tech	6.3744	A
Year Fourth
Year Fourth	Hill Farm	4.6541	B

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Seventh	1	68.68	80.17	<.0001

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Location	_Location	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Seventh	Hill Farm	LA Tech	-1.6756	0.1871	68.68	-8.95	<.0001	<.0001

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Location	Estimate
Year Seventh	LA Tech	12.2550	A
Year Seventh
Year Seventh	Hill Farm	10.5794	B

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Sixth	1	68.68	59.98	<.0001

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Location	_Location	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Sixth	Hill Farm	LA Tech	-1.4494	0.1871	68.68	-7.74	<.0001	<.0001

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Location	Estimate
Year Sixth	LA Tech	11.0402	A
Year Sixth
Year Sixth	Hill Farm	9.5908	B

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Third	1	68.68	12.82	0.0006

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Location	_Location	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Third	Hill Farm	LA Tech	-0.6700	0.1871	68.68	-3.58	0.0006	0.0005

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Location	Estimate
Year Third	LA Tech	2.9246	A
Year Third
Year Third	Hill Farm	2.2546	B

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Fifth	5	68.68	14.17	<.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Family	_Family	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Fifth	AGH2	AGH25	-0.4203	0.3241	68.68	-1.30	0.1991	0.7865
Year Fifth	AGH2	AGHS1	-0.9569	0.3241	68.68	-2.95	0.0043	0.0420
Year Fifth	AGH2	AGHS2	-1.1507	0.3241	68.68	-3.55	0.0007	0.0068
Year Fifth	AGH2	AGHS3	-1.2211	0.3241	68.68	-3.77	0.0003	0.0032
Year Fifth	AGH2	AGHS4	-2.5250	0.3241	68.68	-7.79	<.0001	<.0001
Year Fifth	AGH25	AGHS1	-0.5366	0.3241	68.68	-1.66	0.1024	0.5634
Year Fifth	AGH25	AGHS2	-0.7304	0.3241	68.68	-2.25	0.0274	0.2204
Year Fifth	AGH25	AGHS3	-0.8007	0.3241	68.68	-2.47	0.0160	0.1400
Year Fifth	AGH25	AGHS4	-2.1047	0.3241	68.68	-6.49	<.0001	<.0001
Year Fifth	AGHS1	AGHS2	-0.1938	0.3241	68.68	-0.60	0.5519	0.9910
Year Fifth	AGHS1	AGHS3	-0.2641	0.3241	68.68	-0.81	0.4179	0.9644
Year Fifth	AGHS1	AGHS4	-1.5681	0.3241	68.68	-4.84	<.0001	<.0001
Year Fifth	AGHS2	AGHS3	-0.07037	0.3241	68.68	-0.22	0.8288	0.9999
Year Fifth	AGHS2	AGHS4	-1.3743	0.3241	68.68	-4.24	<.0001	0.0006
Year Fifth	AGHS3	AGHS4	-1.3040	0.3241	68.68	-4.02	0.0001	0.0013

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Family	Estimate
Year Fifth	AGHS4	9.5920		A
Year Fifth
Year Fifth	AGHS3	8.2881		B
Year Fifth				B
Year Fifth	AGHS2	8.2177		B
Year Fifth				B
Year Fifth	AGHS1	8.0239		B
Year Fifth				B
Year Fifth	AGH25	7.4873	C	B
Year Fifth			C
Year Fifth	AGH2	7.0670	C

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Fourth	5	68.68	16.51	<.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Family	_Family	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Fourth	AGH2	AGH25	-0.3639	0.3241	68.68	-1.12	0.2654	0.8712
Year Fourth	AGH2	AGHS1	-1.6194	0.3241	68.68	-5.00	<.0001	<.0001
Year Fourth	AGH2	AGHS2	-1.6211	0.3241	68.68	-5.00	<.0001	<.0001
Year Fourth	AGH2	AGHS3	-1.1847	0.3241	68.68	-3.65	0.0005	0.0048
Year Fourth	AGH2	AGHS4	-2.5590	0.3241	68.68	-7.89	<.0001	<.0001
Year Fourth	AGH25	AGHS1	-1.2555	0.3241	68.68	-3.87	0.0002	0.0022
Year Fourth	AGH25	AGHS2	-1.2572	0.3241	68.68	-3.88	0.0002	0.0022
Year Fourth	AGH25	AGHS3	-0.8208	0.3241	68.68	-2.53	0.0136	0.1217
Year Fourth	AGH25	AGHS4	-2.1950	0.3241	68.68	-6.77	<.0001	<.0001
Year Fourth	AGHS1	AGHS2	-0.00171	0.3241	68.68	-0.01	0.9958	1.0000
Year Fourth	AGHS1	AGHS3	0.4347	0.3241	68.68	1.34	0.1843	0.7617
Year Fourth	AGHS1	AGHS4	-0.9396	0.3241	68.68	-2.90	0.0050	0.0487
Year Fourth	AGHS2	AGHS3	0.4364	0.3241	68.68	1.35	0.1826	0.7586
Year Fourth	AGHS2	AGHS4	-0.9379	0.3241	68.68	-2.89	0.0051	0.0493
Year Fourth	AGHS3	AGHS4	-1.3742	0.3241	68.68	-4.24	<.0001	0.0006

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Family	Estimate
Year Fourth	AGHS4	6.8485		A
Year Fourth
Year Fourth	AGHS2	5.9107		B
Year Fourth				B
Year Fourth	AGHS1	5.9090		B
Year Fourth				B
Year Fourth	AGHS3	5.4743	C	B
Year Fourth			C
Year Fourth	AGH25	4.6535	C	D
Year Fourth				D
Year Fourth	AGH2	4.2896		D

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Seventh	5	68.68	13.26	<.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Family	_Family	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Seventh	AGH2	AGH25	-0.02369	0.3241	68.68	-0.07	0.9420	1.0000
Year Seventh	AGH2	AGHS1	-0.4087	0.3241	68.68	-1.26	0.2116	0.8056
Year Seventh	AGH2	AGHS2	-0.7000	0.3241	68.68	-2.16	0.0343	0.2634
Year Seventh	AGH2	AGHS3	-0.5881	0.3241	68.68	-1.81	0.0740	0.4599
Year Seventh	AGH2	AGHS4	-2.2643	0.3241	68.68	-6.99	<.0001	<.0001
Year Seventh	AGH25	AGHS1	-0.3850	0.3241	68.68	-1.19	0.2390	0.8420
Year Seventh	AGH25	AGHS2	-0.6763	0.3241	68.68	-2.09	0.0407	0.3003
Year Seventh	AGH25	AGHS3	-0.5644	0.3241	68.68	-1.74	0.0861	0.5070
Year Seventh	AGH25	AGHS4	-2.2406	0.3241	68.68	-6.91	<.0001	<.0001
Year Seventh	AGHS1	AGHS2	-0.2912	0.3241	68.68	-0.90	0.3721	0.9463
Year Seventh	AGHS1	AGHS3	-0.1794	0.3241	68.68	-0.55	0.5818	0.9937
Year Seventh	AGHS1	AGHS4	-1.8555	0.3241	68.68	-5.72	<.0001	<.0001
Year Seventh	AGHS2	AGHS3	0.1119	0.3241	68.68	0.35	0.7310	0.9993
Year Seventh	AGHS2	AGHS4	-1.5643	0.3241	68.68	-4.83	<.0001	<.0001
Year Seventh	AGHS3	AGHS4	-1.6762	0.3241	68.68	-5.17	<.0001	<.0001

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Family	Estimate
Year Seventh	AGHS4	13.0174	A
Year Seventh
Year Seventh	AGHS2	11.4531	B
Year Seventh			B
Year Seventh	AGHS3	11.3412	B
Year Seventh			B
Year Seventh	AGHS1	11.1618	B
Year Seventh			B
Year Seventh	AGH25	10.7768	B
Year Seventh			B
Year Seventh	AGH2	10.7531	B

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Sixth	5	68.68	16.04	<.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Family	_Family	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Sixth	AGH2	AGH25	-0.2837	0.3241	68.68	-0.88	0.3845	0.9519
Year Sixth	AGH2	AGHS1	-0.9029	0.3241	68.68	-2.79	0.0069	0.0656
Year Sixth	AGH2	AGHS2	-0.9487	0.3241	68.68	-2.93	0.0046	0.0451
Year Sixth	AGH2	AGHS3	-1.3610	0.3241	68.68	-4.20	<.0001	0.0007
Year Sixth	AGH2	AGHS4	-2.5997	0.3241	68.68	-8.02	<.0001	<.0001
Year Sixth	AGH25	AGHS1	-0.6192	0.3241	68.68	-1.91	0.0603	0.4002
Year Sixth	AGH25	AGHS2	-0.6650	0.3241	68.68	-2.05	0.0440	0.3187
Year Sixth	AGH25	AGHS3	-1.0774	0.3241	68.68	-3.32	0.0014	0.0141
Year Sixth	AGH25	AGHS4	-2.3160	0.3241	68.68	-7.15	<.0001	<.0001
Year Sixth	AGHS1	AGHS2	-0.04586	0.3241	68.68	-0.14	0.8879	1.0000
Year Sixth	AGHS1	AGHS3	-0.4582	0.3241	68.68	-1.41	0.1620	0.7188
Year Sixth	AGHS1	AGHS4	-1.6969	0.3241	68.68	-5.23	<.0001	<.0001
Year Sixth	AGHS2	AGHS3	-0.4123	0.3241	68.68	-1.27	0.2077	0.7998
Year Sixth	AGHS2	AGHS4	-1.6510	0.3241	68.68	-5.09	<.0001	<.0001
Year Sixth	AGHS3	AGHS4	-1.2387	0.3241	68.68	-3.82	0.0003	0.0027

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Family	Estimate
Year Sixth	AGHS4	11.8992		A
Year Sixth
Year Sixth	AGHS3	10.6605		B
Year Sixth				B
Year Sixth	AGHS2	10.2482	C	B
Year Sixth			C	B
Year Sixth	AGHS1	10.2023	C	B	D
Year Sixth			C		D
Year Sixth	AGH25	9.5831	C		D
Year Sixth					D
Year Sixth	AGH2	9.2995			D

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Third	5	68.68	12.50	<.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Family	_Family	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Third	AGH2	AGH25	-0.09482	0.3241	68.68	-0.29	0.7708	0.9997
Year Third	AGH2	AGHS1	-1.5896	0.3241	68.68	-4.90	<.0001	<.0001
Year Third	AGH2	AGHS2	-1.6557	0.3241	68.68	-5.11	<.0001	<.0001
Year Third	AGH2	AGHS3	-0.6881	0.3241	68.68	-2.12	0.0374	0.2815
Year Third	AGH2	AGHS4	-1.7722	0.3241	68.68	-5.47	<.0001	<.0001
Year Third	AGH25	AGHS1	-1.4948	0.3241	68.68	-4.61	<.0001	0.0001
Year Third	AGH25	AGHS2	-1.5609	0.3241	68.68	-4.82	<.0001	<.0001
Year Third	AGH25	AGHS3	-0.5933	0.3241	68.68	-1.83	0.0715	0.4497
Year Third	AGH25	AGHS4	-1.6774	0.3241	68.68	-5.17	<.0001	<.0001
Year Third	AGHS1	AGHS2	-0.06610	0.3241	68.68	-0.20	0.8390	1.0000
Year Third	AGHS1	AGHS3	0.9015	0.3241	68.68	2.78	0.0070	0.0663
Year Third	AGHS1	AGHS4	-0.1826	0.3241	68.68	-0.56	0.5751	0.9932
Year Third	AGHS2	AGHS3	0.9676	0.3241	68.68	2.99	0.0039	0.0384
Year Third	AGHS2	AGHS4	-0.1165	0.3241	68.68	-0.36	0.7204	0.9992
Year Third	AGHS3	AGHS4	-1.0841	0.3241	68.68	-3.34	0.0013	0.0132

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Family	Estimate
Year Third	AGHS4	3.3951		A
Year Third				A
Year Third	AGHS2	3.2786		A
Year Third				A
Year Third	AGHS1	3.2125	B	A
Year Third			B
Year Third	AGHS3	2.3110	B	C
Year Third				C
Year Third	AGH25	1.7177		C
Year Third				C
Year Third	AGH2	1.6229		C

LS-Means Results

Family LS-Means

Family Least Squares Means
Family	Estimate	Standard Error	DF	t Value	Pr > |t|
AGH2	6.6064	0.1930	36	34.23	<.0001
AGH25	6.8437	0.1930	36	35.45	<.0001
AGHS1	7.7019	0.1930	36	39.90	<.0001
AGHS2	7.8216	0.1930	36	40.52	<.0001
AGHS3	7.6150	0.1930	36	39.45	<.0001
AGHS4	8.9504	0.1930	36	46.37	<.0001

Family Diffs

Differences of Family Least Squares Means Adjustment for Multiple Comparisons: Tukey
Family	_Family	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
AGH2	AGH25	-0.2373	0.2730	36	-0.87	0.3905	0.9514
AGH2	AGHS1	-1.0955	0.2730	36	-4.01	0.0003	0.0037
AGH2	AGHS2	-1.2152	0.2730	36	-4.45	<.0001	0.0010
AGH2	AGHS3	-1.0086	0.2730	36	-3.69	0.0007	0.0088
AGH2	AGHS4	-2.3440	0.2730	36	-8.59	<.0001	<.0001
AGH25	AGHS1	-0.8582	0.2730	36	-3.14	0.0033	0.0362
AGH25	AGHS2	-0.9780	0.2730	36	-3.58	0.0010	0.0118
AGH25	AGHS3	-0.7713	0.2730	36	-2.83	0.0077	0.0760
AGH25	AGHS4	-2.1067	0.2730	36	-7.72	<.0001	<.0001
AGHS1	AGHS2	-0.1197	0.2730	36	-0.44	0.6635	0.9978
AGHS1	AGHS3	0.08690	0.2730	36	0.32	0.7521	0.9995
AGHS1	AGHS4	-1.2485	0.2730	36	-4.57	<.0001	0.0007
AGHS2	AGHS3	0.2066	0.2730	36	0.76	0.4540	0.9729
AGHS2	AGHS4	-1.1288	0.2730	36	-4.14	0.0002	0.0026
AGHS3	AGHS4	-1.3354	0.2730	36	-4.89	<.0001	0.0003

Family Lines

Tukey Grouping for Family Least Squares Means (Alpha=0.05)
LS-means with the same letter are not significantly different.
Family	Estimate
AGHS4	8.9504		A
AGHS2	7.8216		B
			B
AGHS1	7.7019		B
			B
AGHS3	7.6150	C	B
		C
AGH25	6.8437	C	D
			D
AGH2	6.6064		D

SAS Output Height

htmltools::includeHTML("https://raw.githubusercontent.com/azd169/SAS_outputs/main/hane-results-height.html")

## Warning: `includeHTML()` was provided a `path` that appears to be a complete HTML document.
## ✖ Path: https://raw.githubusercontent.com/azd169/SAS_outputs/main/hane-results-height.html
## ℹ Use `tags$iframe()` to include an HTML document. You can either ensure `path` is accessible in your app or document (see e.g. `shiny::addResourcePath()`) and pass the relative path to the `src` argument. Or you can read the contents of `path` and pass the contents to `srcdoc`.

Results: hane.sas

The GLIMMIX Procedure

The Glimmix Procedure

Model Information

Model Information
Data Set	WORK.HEIGHT
Response Variable	Height
Response Distribution	Gaussian
Link Function	Identity
Variance Function	Default
Variance Matrix Blocked By	GLB
Estimation Technique	Restricted Maximum Likelihood
Degrees of Freedom Method	Kenward-Roger
Fixed Effects SE Adjustment	Kenward-Roger

Class Level Information

Class Level Information
Class	Levels	Values
Year	8	Fifth First Fourth Initial Second Seventh Sixth Third
Location	2	Hill Farm LA Tech
Block	4	1 2 3 4
Family	6	AGH2 AGH25 AGHS1 AGHS2 AGHS3 AGHS4
GLB	48	AGH25_Hill Farm_1 AGH25_Hill Farm_2 AGH25_Hill Farm_3 AGH25_Hill Farm_4 AGH25_LA Tech_1 AGH25_LA Tech_2 AGH25_LA Tech_3 AGH25_LA Tech_4 AGH2_Hill Farm_1 AGH2_Hill Farm_2 AGH2_Hill Farm_3 AGH2_Hill Farm_4 AGH2_LA Tech_1 AGH2_LA Tech_2 AGH2_LA Tech_3 AGH2_LA Tech_4 AGHS1_Hill Farm_1 AGHS1_Hill Farm_2 AGHS1_Hill Farm_3 AGHS1_Hill Farm_4 AGHS1_LA Tech_1 AGHS1_LA Tech_2 AGHS1_LA Tech_3 AGHS1_LA Tech_4 AGHS2_Hill Farm_1 AGHS2_Hill Farm_2 AGHS2_Hill Farm_3 AGHS2_Hill Farm_4 AGHS2_LA Tech_1 AGHS2_LA Tech_2 AGHS2_LA Tech_3 AGHS2_LA Tech_4 AGHS3_Hill Farm_1 AGHS3_Hill Farm_2 AGHS3_Hill Farm_3 AGHS3_Hill Farm_4 AGHS3_LA Tech_1 AGHS3_LA Tech_2 AGHS3_LA Tech_3 AGHS3_LA Tech_4 AGHS4_Hill Farm_1 AGHS4_Hill Farm_2 AGHS4_Hill Farm_3 AGHS4_Hill Farm_4 AGHS4_LA Tech_1 AGHS4_LA Tech_2 AGHS4_LA Tech_3 AGHS4_LA Tech_4

Number of Observations

Number of Observations Read	384
Number of Observations Used	384

Dimensions

Dimensions
G-side Cov. Parameters	2
R-side Cov. Parameters	1
Columns in X	189
Columns in Z per Subject	1
Subjects (Blocks in V)	48
Max Obs per Subject	8

Optimization Information

Optimization Information
Optimization Technique	Dual Quasi-Newton
Parameters in Optimization	2
Lower Boundaries	2
Upper Boundaries	1
Fixed Effects	Profiled
Residual Variance	Profiled
Starting From	Data

Iteration History

Iteration History
Iteration	Restarts	Evaluations	Objective Function	Change	Max Gradient
0	0	4	419.85607123	.	3.56E-12

Convergence Status

Convergence criterion (ABSGCONV=0.00001) satisfied.

Hessian Matrix

Hessian Matrix
Cov Parm	Subject	Row	Col1	Col2
Variance	GLB	1	80.0040
AR(1)	GLB	2

Fit Statistics

Fit Statistics
-2 Res Log Likelihood	419.86
AIC (smaller is better)	425.86
AICC (smaller is better)	425.94
BIC (smaller is better)	431.47
CAIC (smaller is better)	434.47
HQIC (smaller is better)	427.98
Generalized Chi-Square	37.31
Gener. Chi-Square / DF	0.13

Covariance Parameter Estimates

Covariance Parameter Estimates
Cov Parm	Subject	Estimate	Standard Error
Variance	GLB	0.06509	0.01921
AR(1)	GLB	0	.
Residual		0.1295	0.01154

Solutions for Fixed Effects

Solutions for Fixed Effects
Effect	Year	Location	Family	Estimate	Standard Error	DF	t Value	Pr > |t|
Intercept				3.6033	0.2206	161.5	16.34	<.0001
Year	Fifth			3.7856	0.2545	252	14.88	<.0001
Year	First			-2.6121	0.2545	252	-10.26	<.0001
Year	Fourth			2.2226	0.2545	252	8.73	<.0001
Year	Initial			-3.0528	0.2545	252	-12.00	<.0001
Year	Second			-1.2216	0.2545	252	-4.80	<.0001
Year	Seventh			5.8582	0.2545	252	23.02	<.0001
Year	Sixth			5.9387	0.2545	252	23.34	<.0001
Year	Third			0	.	.	.	.
Location		Hill Farm		-0.7677	0.3119	161.5	-2.46	0.0149
Location		LA Tech		0	.	.	.	.
Year*Location	Fifth	Hill Farm		-0.3317	0.3599	252	-0.92	0.3575
Year*Location	Fifth	LA Tech		0	.	.	.	.
Year*Location	First	Hill Farm		0.7031	0.3599	252	1.95	0.0518
Year*Location	First	LA Tech		0	.	.	.	.
Year*Location	Fourth	Hill Farm		-0.5411	0.3599	252	-1.50	0.1340
Year*Location	Fourth	LA Tech		0	.	.	.	.
Year*Location	Initial	Hill Farm		0.7707	0.3599	252	2.14	0.0332
Year*Location	Initial	LA Tech		0	.	.	.	.
Year*Location	Second	Hill Farm		0.4508	0.3599	252	1.25	0.2115
Year*Location	Second	LA Tech		0	.	.	.	.
Year*Location	Seventh	Hill Farm		1.0020	0.3599	252	2.78	0.0058
Year*Location	Seventh	LA Tech		0	.	.	.	.
Year*Location	Sixth	Hill Farm		-2.5706	0.3599	252	-7.14	<.0001
Year*Location	Sixth	LA Tech		0	.	.	.	.
Year*Location	Third	Hill Farm		0	.	.	.	.
Year*Location	Third	LA Tech		0	.	.	.	.
Family			AGH2	-0.9547	0.3119	161.5	-3.06	0.0026
Family			AGH25	-0.7602	0.3119	161.5	-2.44	0.0159
Family			AGHS1	0.2949	0.3119	161.5	0.95	0.3459
Family			AGHS2	0.4191	0.3119	161.5	1.34	0.1810
Family			AGHS3	-0.3812	0.3119	161.5	-1.22	0.2235
Family			AGHS4	0	.	.	.	.
Year*Family	Fifth		AGH2	-0.4696	0.3599	252	-1.30	0.1931
Year*Family	Fifth		AGH25	-0.04786	0.3599	252	-0.13	0.8943
Year*Family	Fifth		AGHS1	-0.1533	0.3599	252	-0.43	0.6706
Year*Family	Fifth		AGHS2	-0.05521	0.3599	252	-0.15	0.8782
Year*Family	Fifth		AGHS3	-0.09546	0.3599	252	-0.27	0.7910
Year*Family	Fifth		AGHS4	0	.	.	.	.
Year*Family	First		AGH2	0.6650	0.3599	252	1.85	0.0658
Year*Family	First		AGH25	0.4395	0.3599	252	1.22	0.2232
Year*Family	First		AGHS1	-0.3370	0.3599	252	-0.94	0.3500
Year*Family	First		AGHS2	-0.4107	0.3599	252	-1.14	0.2549
Year*Family	First		AGHS3	0.2685	0.3599	252	0.75	0.4563
Year*Family	First		AGHS4	0	.	.	.	.
Year*Family	Fourth		AGH2	-0.5640	0.3599	252	-1.57	0.1184
Year*Family	Fourth		AGH25	-0.3129	0.3599	252	-0.87	0.3855
Year*Family	Fourth		AGHS1	-0.06452	0.3599	252	-0.18	0.8579
Year*Family	Fourth		AGHS2	0.07452	0.3599	252	0.21	0.8361
Year*Family	Fourth		AGHS3	0.06692	0.3599	252	0.19	0.8526
Year*Family	Fourth		AGHS4	0	.	.	.	.
Year*Family	Initial		AGH2	0.8553	0.3599	252	2.38	0.0182
Year*Family	Initial		AGH25	0.6140	0.3599	252	1.71	0.0892
Year*Family	Initial		AGHS1	-0.3388	0.3599	252	-0.94	0.3474
Year*Family	Initial		AGHS2	-0.5205	0.3599	252	-1.45	0.1494
Year*Family	Initial		AGHS3	0.3439	0.3599	252	0.96	0.3402
Year*Family	Initial		AGHS4	0	.	.	.	.
Year*Family	Second		AGH2	0.1622	0.3599	252	0.45	0.6525
Year*Family	Second		AGH25	0.1249	0.3599	252	0.35	0.7288
Year*Family	Second		AGHS1	-0.1838	0.3599	252	-0.51	0.6100
Year*Family	Second		AGHS2	-0.05966	0.3599	252	-0.17	0.8685
Year*Family	Second		AGHS3	0.1786	0.3599	252	0.50	0.6201
Year*Family	Second		AGHS4	0	.	.	.	.
Year*Family	Seventh		AGH2	-0.1003	0.3599	252	-0.28	0.7808
Year*Family	Seventh		AGH25	-0.1628	0.3599	252	-0.45	0.6514
Year*Family	Seventh		AGHS1	-0.2814	0.3599	252	-0.78	0.4351
Year*Family	Seventh		AGHS2	-0.3133	0.3599	252	-0.87	0.3849
Year*Family	Seventh		AGHS3	0.3029	0.3599	252	0.84	0.4009
Year*Family	Seventh		AGHS4	0	.	.	.	.
Year*Family	Sixth		AGH2	-1.4486	0.3599	252	-4.02	<.0001
Year*Family	Sixth		AGH25	-0.9490	0.3599	252	-2.64	0.0089
Year*Family	Sixth		AGHS1	-0.06103	0.3599	252	-0.17	0.8655
Year*Family	Sixth		AGHS2	-0.1143	0.3599	252	-0.32	0.7511
Year*Family	Sixth		AGHS3	-0.05486	0.3599	252	-0.15	0.8790
Year*Family	Sixth		AGHS4	0	.	.	.	.
Year*Family	Third		AGH2	0	.	.	.	.
Year*Family	Third		AGH25	0	.	.	.	.
Year*Family	Third		AGHS1	0	.	.	.	.
Year*Family	Third		AGHS2	0	.	.	.	.
Year*Family	Third		AGHS3	0	.	.	.	.
Year*Family	Third		AGHS4	0	.	.	.	.
Location*Family		Hill Farm	AGH2	0.4218	0.4412	161.5	0.96	0.3404
Location*Family		Hill Farm	AGH25	0.3009	0.4412	161.5	0.68	0.4962
Location*Family		Hill Farm	AGHS1	-0.02870	0.4412	161.5	-0.07	0.9482
Location*Family		Hill Farm	AGHS2	-0.3231	0.4412	161.5	-0.73	0.4650
Location*Family		Hill Farm	AGHS3	-0.05930	0.4412	161.5	-0.13	0.8932
Location*Family		Hill Farm	AGHS4	0	.	.	.	.
Location*Family		LA Tech	AGH2	0	.	.	.	.
Location*Family		LA Tech	AGH25	0	.	.	.	.
Location*Family		LA Tech	AGHS1	0	.	.	.	.
Location*Family		LA Tech	AGHS2	0	.	.	.	.
Location*Family		LA Tech	AGHS3	0	.	.	.	.
Location*Family		LA Tech	AGHS4	0	.	.	.	.
Year*Location*Family	Fifth	Hill Farm	AGH2	-0.4498	0.5090	252	-0.88	0.3776
Year*Location*Family	Fifth	Hill Farm	AGH25	-0.5438	0.5090	252	-1.07	0.2863
Year*Location*Family	Fifth	Hill Farm	AGHS1	-0.3138	0.5090	252	-0.62	0.5381
Year*Location*Family	Fifth	Hill Farm	AGHS2	0.1410	0.5090	252	0.28	0.7820
Year*Location*Family	Fifth	Hill Farm	AGHS3	-0.3825	0.5090	252	-0.75	0.4531
Year*Location*Family	Fifth	Hill Farm	AGHS4	0	.	.	.	.
Year*Location*Family	Fifth	LA Tech	AGH2	0	.	.	.	.
Year*Location*Family	Fifth	LA Tech	AGH25	0	.	.	.	.
Year*Location*Family	Fifth	LA Tech	AGHS1	0	.	.	.	.
Year*Location*Family	Fifth	LA Tech	AGHS2	0	.	.	.	.
Year*Location*Family	Fifth	LA Tech	AGHS3	0	.	.	.	.
Year*Location*Family	Fifth	LA Tech	AGHS4	0	.	.	.	.
Year*Location*Family	First	Hill Farm	AGH2	-0.3662	0.5090	252	-0.72	0.4725
Year*Location*Family	First	Hill Farm	AGH25	-0.1775	0.5090	252	-0.35	0.7276
Year*Location*Family	First	Hill Farm	AGHS1	0.1585	0.5090	252	0.31	0.7557
Year*Location*Family	First	Hill Farm	AGHS2	0.3192	0.5090	252	0.63	0.5312
Year*Location*Family	First	Hill Farm	AGHS3	0.1219	0.5090	252	0.24	0.8109
Year*Location*Family	First	Hill Farm	AGHS4	0	.	.	.	.
Year*Location*Family	First	LA Tech	AGH2	0	.	.	.	.
Year*Location*Family	First	LA Tech	AGH25	0	.	.	.	.
Year*Location*Family	First	LA Tech	AGHS1	0	.	.	.	.
Year*Location*Family	First	LA Tech	AGHS2	0	.	.	.	.
Year*Location*Family	First	LA Tech	AGHS3	0	.	.	.	.
Year*Location*Family	First	LA Tech	AGHS4	0	.	.	.	.
Year*Location*Family	Fourth	Hill Farm	AGH2	-0.1764	0.5090	252	-0.35	0.7292
Year*Location*Family	Fourth	Hill Farm	AGH25	-0.2610	0.5090	252	-0.51	0.6086
Year*Location*Family	Fourth	Hill Farm	AGHS1	-0.08584	0.5090	252	-0.17	0.8662
Year*Location*Family	Fourth	Hill Farm	AGHS2	0.07420	0.5090	252	0.15	0.8842
Year*Location*Family	Fourth	Hill Farm	AGHS3	-0.2134	0.5090	252	-0.42	0.6754
Year*Location*Family	Fourth	Hill Farm	AGHS4	0	.	.	.	.
Year*Location*Family	Fourth	LA Tech	AGH2	0	.	.	.	.
Year*Location*Family	Fourth	LA Tech	AGH25	0	.	.	.	.
Year*Location*Family	Fourth	LA Tech	AGHS1	0	.	.	.	.
Year*Location*Family	Fourth	LA Tech	AGHS2	0	.	.	.	.
Year*Location*Family	Fourth	LA Tech	AGHS3	0	.	.	.	.
Year*Location*Family	Fourth	LA Tech	AGHS4	0	.	.	.	.
Year*Location*Family	Initial	Hill Farm	AGH2	-0.4515	0.5090	252	-0.89	0.3759
Year*Location*Family	Initial	Hill Farm	AGH25	-0.2910	0.5090	252	-0.57	0.5681
Year*Location*Family	Initial	Hill Farm	AGHS1	0.08518	0.5090	252	0.17	0.8672
Year*Location*Family	Initial	Hill Farm	AGHS2	0.3375	0.5090	252	0.66	0.5079
Year*Location*Family	Initial	Hill Farm	AGHS3	0.08876	0.5090	252	0.17	0.8617
Year*Location*Family	Initial	Hill Farm	AGHS4	0	.	.	.	.
Year*Location*Family	Initial	LA Tech	AGH2	0	.	.	.	.
Year*Location*Family	Initial	LA Tech	AGH25	0	.	.	.	.
Year*Location*Family	Initial	LA Tech	AGHS1	0	.	.	.	.
Year*Location*Family	Initial	LA Tech	AGHS2	0	.	.	.	.
Year*Location*Family	Initial	LA Tech	AGHS3	0	.	.	.	.
Year*Location*Family	Initial	LA Tech	AGHS4	0	.	.	.	.
Year*Location*Family	Second	Hill Farm	AGH2	-0.1221	0.5090	252	-0.24	0.8106
Year*Location*Family	Second	Hill Farm	AGH25	-0.1298	0.5090	252	-0.26	0.7989
Year*Location*Family	Second	Hill Farm	AGHS1	0.06376	0.5090	252	0.13	0.9004
Year*Location*Family	Second	Hill Farm	AGHS2	0.1648	0.5090	252	0.32	0.7463
Year*Location*Family	Second	Hill Farm	AGHS3	0.1828	0.5090	252	0.36	0.7197
Year*Location*Family	Second	Hill Farm	AGHS4	0	.	.	.	.
Year*Location*Family	Second	LA Tech	AGH2	0	.	.	.	.
Year*Location*Family	Second	LA Tech	AGH25	0	.	.	.	.
Year*Location*Family	Second	LA Tech	AGHS1	0	.	.	.	.
Year*Location*Family	Second	LA Tech	AGHS2	0	.	.	.	.
Year*Location*Family	Second	LA Tech	AGHS3	0	.	.	.	.
Year*Location*Family	Second	LA Tech	AGHS4	0	.	.	.	.
Year*Location*Family	Seventh	Hill Farm	AGH2	-0.6213	0.5090	252	-1.22	0.2234
Year*Location*Family	Seventh	Hill Farm	AGH25	-0.9667	0.5090	252	-1.90	0.0587
Year*Location*Family	Seventh	Hill Farm	AGHS1	0.1514	0.5090	252	0.30	0.7664
Year*Location*Family	Seventh	Hill Farm	AGHS2	0.2505	0.5090	252	0.49	0.6230
Year*Location*Family	Seventh	Hill Farm	AGHS3	-1.0173	0.5090	252	-2.00	0.0467
Year*Location*Family	Seventh	Hill Farm	AGHS4	0	.	.	.	.
Year*Location*Family	Seventh	LA Tech	AGH2	0	.	.	.	.
Year*Location*Family	Seventh	LA Tech	AGH25	0	.	.	.	.
Year*Location*Family	Seventh	LA Tech	AGHS1	0	.	.	.	.
Year*Location*Family	Seventh	LA Tech	AGHS2	0	.	.	.	.
Year*Location*Family	Seventh	LA Tech	AGHS3	0	.	.	.	.
Year*Location*Family	Seventh	LA Tech	AGHS4	0	.	.	.	.
Year*Location*Family	Sixth	Hill Farm	AGH2	1.1240	0.5090	252	2.21	0.0281
Year*Location*Family	Sixth	Hill Farm	AGH25	0.5276	0.5090	252	1.04	0.3010
Year*Location*Family	Sixth	Hill Farm	AGHS1	0.06816	0.5090	252	0.13	0.8936
Year*Location*Family	Sixth	Hill Farm	AGHS2	0.3570	0.5090	252	0.70	0.4838
Year*Location*Family	Sixth	Hill Farm	AGHS3	0.1141	0.5090	252	0.22	0.8228
Year*Location*Family	Sixth	Hill Farm	AGHS4	0	.	.	.	.
Year*Location*Family	Sixth	LA Tech	AGH2	0	.	.	.	.
Year*Location*Family	Sixth	LA Tech	AGH25	0	.	.	.	.
Year*Location*Family	Sixth	LA Tech	AGHS1	0	.	.	.	.
Year*Location*Family	Sixth	LA Tech	AGHS2	0	.	.	.	.
Year*Location*Family	Sixth	LA Tech	AGHS3	0	.	.	.	.
Year*Location*Family	Sixth	LA Tech	AGHS4	0	.	.	.	.
Year*Location*Family	Third	Hill Farm	AGH2	0	.	.	.	.
Year*Location*Family	Third	Hill Farm	AGH25	0	.	.	.	.
Year*Location*Family	Third	Hill Farm	AGHS1	0	.	.	.	.
Year*Location*Family	Third	Hill Farm	AGHS2	0	.	.	.	.
Year*Location*Family	Third	Hill Farm	AGHS3	0	.	.	.	.
Year*Location*Family	Third	Hill Farm	AGHS4	0	.	.	.	.
Year*Location*Family	Third	LA Tech	AGH2	0	.	.	.	.
Year*Location*Family	Third	LA Tech	AGH25	0	.	.	.	.
Year*Location*Family	Third	LA Tech	AGHS1	0	.	.	.	.
Year*Location*Family	Third	LA Tech	AGHS2	0	.	.	.	.
Year*Location*Family	Third	LA Tech	AGHS3	0	.	.	.	.
Year*Location*Family	Third	LA Tech	AGHS4	0	.	.	.	.

Type III Tests of Fixed Effects

Type III Tests of Fixed Effects
Effect	Num DF	Den DF	F Value	Pr > F
Year	7	252	3715.43	<.0001
Location	1	36	100.98	<.0001
Year*Location	7	252	95.45	<.0001
Family	5	36	21.61	<.0001
Year*Family	35	252	5.27	<.0001
Location*Family	5	36	0.70	0.6304
Year*Location*Family	35	252	1.24	0.1731

Solution for Random Effects

Solution for Random Effects
Effect	Subject	Estimate	Std Err Pred	DF	t Value	Pr > |t|	Alpha	Lower	Upper
Intercept	GLB AGH25_Hill Farm_1	-0.1321	0.1622	46.92	-0.81	0.4196	0.05	-0.4584	0.1942
Intercept	GLB AGH25_Hill Farm_2	0.06641	0.1622	46.92	0.41	0.6841	0.05	-0.2599	0.3927
Intercept	GLB AGH25_Hill Farm_3	-0.09112	0.1622	46.92	-0.56	0.5769	0.05	-0.4174	0.2352
Intercept	GLB AGH25_Hill Farm_4	0.1568	0.1622	46.92	0.97	0.3387	0.05	-0.1695	0.4831
Intercept	GLB AGH25_LA Tech_1	-0.1835	0.1622	46.92	-1.13	0.2637	0.05	-0.5098	0.1428
Intercept	GLB AGH25_LA Tech_2	-0.1597	0.1622	46.92	-0.98	0.3299	0.05	-0.4860	0.1666
Intercept	GLB AGH25_LA Tech_3	0.07326	0.1622	46.92	0.45	0.6536	0.05	-0.2530	0.3995
Intercept	GLB AGH25_LA Tech_4	0.2699	0.1622	46.92	1.66	0.1027	0.05	-0.05639	0.5962
Intercept	GLB AGH2_Hill Farm_1	-0.07375	0.1622	46.92	-0.45	0.6514	0.05	-0.4000	0.2525
Intercept	GLB AGH2_Hill Farm_2	-0.08202	0.1622	46.92	-0.51	0.6154	0.05	-0.4083	0.2443
Intercept	GLB AGH2_Hill Farm_3	0.1422	0.1622	46.92	0.88	0.3849	0.05	-0.1841	0.4685
Intercept	GLB AGH2_Hill Farm_4	0.01353	0.1622	46.92	0.08	0.9339	0.05	-0.3128	0.3398
Intercept	GLB AGH2_LA Tech_1	-0.1325	0.1622	46.92	-0.82	0.4180	0.05	-0.4588	0.1938
Intercept	GLB AGH2_LA Tech_2	-0.1581	0.1622	46.92	-0.97	0.3346	0.05	-0.4844	0.1682
Intercept	GLB AGH2_LA Tech_3	0.4081	0.1622	46.92	2.52	0.0153	0.05	0.08183	0.7344
Intercept	GLB AGH2_LA Tech_4	-0.1175	0.1622	46.92	-0.72	0.4724	0.05	-0.4438	0.2088
Intercept	GLB AGHS1_Hill Farm_1	0.03870	0.1622	46.92	0.24	0.8124	0.05	-0.2876	0.3650
Intercept	GLB AGHS1_Hill Farm_2	0.03476	0.1622	46.92	0.21	0.8312	0.05	-0.2915	0.3611
Intercept	GLB AGHS1_Hill Farm_3	0.2445	0.1622	46.92	1.51	0.1384	0.05	-0.08179	0.5708
Intercept	GLB AGHS1_Hill Farm_4	-0.3180	0.1622	46.92	-1.96	0.0559	0.05	-0.6443	0.008324
Intercept	GLB AGHS1_LA Tech_1	-0.2625	0.1622	46.92	-1.62	0.1123	0.05	-0.5888	0.06382
Intercept	GLB AGHS1_LA Tech_2	0.03883	0.1622	46.92	0.24	0.8118	0.05	-0.2875	0.3651
Intercept	GLB AGHS1_LA Tech_3	0.3097	0.1622	46.92	1.91	0.0623	0.05	-0.01662	0.6360
Intercept	GLB AGHS1_LA Tech_4	-0.08603	0.1622	46.92	-0.53	0.5983	0.05	-0.4123	0.2403
Intercept	GLB AGHS2_Hill Farm_1	-0.2649	0.1622	46.92	-1.63	0.1091	0.05	-0.5912	0.06140
Intercept	GLB AGHS2_Hill Farm_2	-0.2065	0.1622	46.92	-1.27	0.2093	0.05	-0.5328	0.1198
Intercept	GLB AGHS2_Hill Farm_3	0.2849	0.1622	46.92	1.76	0.0856	0.05	-0.04144	0.6111
Intercept	GLB AGHS2_Hill Farm_4	0.1865	0.1622	46.92	1.15	0.2560	0.05	-0.1398	0.5128
Intercept	GLB AGHS2_LA Tech_1	-0.2952	0.1622	46.92	-1.82	0.0751	0.05	-0.6215	0.03107
Intercept	GLB AGHS2_LA Tech_2	-0.01269	0.1622	46.92	-0.08	0.9380	0.05	-0.3390	0.3136
Intercept	GLB AGHS2_LA Tech_3	0.3123	0.1622	46.92	1.93	0.0602	0.05	-0.01401	0.6386
Intercept	GLB AGHS2_LA Tech_4	-0.00437	0.1622	46.92	-0.03	0.9786	0.05	-0.3307	0.3219
Intercept	GLB AGHS3_Hill Farm_1	-0.00285	0.1622	46.92	-0.02	0.9861	0.05	-0.3291	0.3234
Intercept	GLB AGHS3_Hill Farm_2	-0.4217	0.1622	46.92	-2.60	0.0124	0.05	-0.7480	-0.09539
Intercept	GLB AGHS3_Hill Farm_3	0.3871	0.1622	46.92	2.39	0.0211	0.05	0.06081	0.7134
Intercept	GLB AGHS3_Hill Farm_4	0.03743	0.1622	46.92	0.23	0.8185	0.05	-0.2889	0.3637
Intercept	GLB AGHS3_LA Tech_1	-0.2095	0.1622	46.92	-1.29	0.2029	0.05	-0.5357	0.1168
Intercept	GLB AGHS3_LA Tech_2	-0.09466	0.1622	46.92	-0.58	0.5623	0.05	-0.4209	0.2316
Intercept	GLB AGHS3_LA Tech_3	0.1860	0.1622	46.92	1.15	0.2573	0.05	-0.1403	0.5123
Intercept	GLB AGHS3_LA Tech_4	0.1181	0.1622	46.92	0.73	0.4701	0.05	-0.2082	0.4444
Intercept	GLB AGHS4_Hill Farm_1	-0.2830	0.1622	46.92	-1.75	0.0875	0.05	-0.6093	0.04324
Intercept	GLB AGHS4_Hill Farm_2	0.1853	0.1622	46.92	1.14	0.2590	0.05	-0.1410	0.5116
Intercept	GLB AGHS4_Hill Farm_3	-0.00433	0.1622	46.92	-0.03	0.9788	0.05	-0.3306	0.3220
Intercept	GLB AGHS4_Hill Farm_4	0.1021	0.1622	46.92	0.63	0.5322	0.05	-0.2242	0.4284
Intercept	GLB AGHS4_LA Tech_1	-0.2516	0.1622	46.92	-1.55	0.1275	0.05	-0.5779	0.07465
Intercept	GLB AGHS4_LA Tech_2	-0.07550	0.1622	46.92	-0.47	0.6437	0.05	-0.4018	0.2508
Intercept	GLB AGHS4_LA Tech_3	0.1680	0.1622	46.92	1.04	0.3055	0.05	-0.1583	0.4943
Intercept	GLB AGHS4_LA Tech_4	0.1591	0.1622	46.92	0.98	0.3316	0.05	-0.1672	0.4854

Slices

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Fifth	1	161.5	105.11	<.0001

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Location	_Location	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Fifth	Hill Farm	LA Tech	-1.3057	0.1274	161.5	-10.25	<.0001	<.0001

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Location	Estimate
Year Fifth	LA Tech	7.0217	A
Year Fifth
Year Fifth	Hill Farm	5.7160	B

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year First	1	161.5	0.00	0.9793

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Location	_Location	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year First	Hill Farm	LA Tech	-0.00330	0.1274	161.5	-0.03	0.9793	0.9793

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Location	Estimate
Year First	LA Tech	0.8650	A
Year First			A
Year First	Hill Farm	0.8617	A

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Fourth	1	161.5	115.26	<.0001

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Location	_Location	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Fourth	Hill Farm	LA Tech	-1.3672	0.1274	161.5	-10.74	<.0001	<.0001

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Location	Estimate
Year Fourth	LA Tech	5.4623	A
Year Fourth
Year Fourth	Hill Farm	4.0950	B

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Initial	1	161.5	0.02	0.8976

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Location	_Location	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Initial	Hill Farm	LA Tech	0.01642	0.1274	161.5	0.13	0.8976	0.8975

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Location	Estimate
Year Initial	Hill Farm	0.4955	A
Year Initial			A
Year Initial	LA Tech	0.4791	A

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Second	1	161.5	3.50	0.0631

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Location	_Location	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Second	Hill Farm	LA Tech	-0.2384	0.1274	161.5	-1.87	0.0631	0.0624

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Location	Estimate
Year Second	LA Tech	2.1884	A
Year Second			A
Year Second	Hill Farm	1.9500	A

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Seventh	1	161.5	0.40	0.5259

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Location	_Location	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Seventh	Hill Farm	LA Tech	-0.08095	0.1274	161.5	-0.64	0.5259	0.5256

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Location	Estimate
Year Seventh	LA Tech	9.1387	A
Year Seventh			A
Year Seventh	Hill Farm	9.0577	A

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Sixth	1	161.5	526.16	<.0001

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Location	_Location	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Sixth	Hill Farm	LA Tech	-2.9212	0.1274	161.5	-22.94	<.0001	<.0001

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Location	Estimate
Year Sixth	LA Tech	8.8736	A
Year Sixth
Year Sixth	Hill Farm	5.9524	B

Year*Location F Test

F Test for Year*Location Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Third	1	161.5	31.59	<.0001

Year*Location Diffs

Simple Differences of Year*Location Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Location	_Location	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Third	Hill Farm	LA Tech	-0.7158	0.1274	161.5	-5.62	<.0001	<.0001

Year*Location Lines

Tukey Grouping for Year*Location Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Location	Estimate
Year Third	LA Tech	3.3729	A
Year Third
Year Third	Hill Farm	2.6572	B

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Fifth	5	161.5	17.82	<.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Family	_Family	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Fifth	AGH2	AGH25	-0.5088	0.2206	161.5	-2.31	0.0223	0.1952
Year Fifth	AGH2	AGHS1	-1.4087	0.2206	161.5	-6.39	<.0001	<.0001
Year Fifth	AGH2	AGHS2	-1.7111	0.2206	161.5	-7.76	<.0001	<.0001
Year Fifth	AGH2	AGHS3	-0.7408	0.2206	161.5	-3.36	0.0010	0.0115
Year Fifth	AGH2	AGHS4	-1.4383	0.2206	161.5	-6.52	<.0001	<.0001
Year Fifth	AGH25	AGHS1	-0.8999	0.2206	161.5	-4.08	<.0001	0.0009
Year Fifth	AGH25	AGHS2	-1.2023	0.2206	161.5	-5.45	<.0001	<.0001
Year Fifth	AGH25	AGHS3	-0.2320	0.2206	161.5	-1.05	0.2945	0.8997
Year Fifth	AGH25	AGHS4	-0.9295	0.2206	161.5	-4.21	<.0001	0.0005
Year Fifth	AGHS1	AGHS2	-0.3024	0.2206	161.5	-1.37	0.1723	0.7443
Year Fifth	AGHS1	AGHS3	0.6679	0.2206	161.5	3.03	0.0029	0.0321
Year Fifth	AGHS1	AGHS4	-0.02959	0.2206	161.5	-0.13	0.8934	1.0000
Year Fifth	AGHS2	AGHS3	0.9704	0.2206	161.5	4.40	<.0001	0.0002
Year Fifth	AGHS2	AGHS4	0.2728	0.2206	161.5	1.24	0.2179	0.8183
Year Fifth	AGHS3	AGHS4	-0.6975	0.2206	161.5	-3.16	0.0019	0.0215

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Family	Estimate
Year Fifth	AGHS2	7.1120		A
Year Fifth				A
Year Fifth	AGHS4	6.8392		A
Year Fifth				A
Year Fifth	AGHS1	6.8096		A
Year Fifth
Year Fifth	AGHS3	6.1416		B
Year Fifth				B
Year Fifth	AGH25	5.9097	C	B
Year Fifth			C
Year Fifth	AGH2	5.4008	C

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year First	5	161.5	0.72	0.6060

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Family	_Family	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year First	AGH2	AGH25	-0.00289	0.2206	161.5	-0.01	0.9896	1.0000
Year First	AGH2	AGHS1	-0.2847	0.2206	161.5	-1.29	0.1986	0.7899
Year First	AGH2	AGHS2	-0.2683	0.2206	161.5	-1.22	0.2256	0.8286
Year First	AGH2	AGHS3	-0.1805	0.2206	161.5	-0.82	0.4143	0.9639
Year First	AGH2	AGHS4	-0.2619	0.2206	161.5	-1.19	0.2368	0.8426
Year First	AGH25	AGHS1	-0.2818	0.2206	161.5	-1.28	0.2032	0.7970
Year First	AGH25	AGHS2	-0.2654	0.2206	161.5	-1.20	0.2306	0.8350
Year First	AGH25	AGHS3	-0.1777	0.2206	161.5	-0.81	0.4218	0.9664
Year First	AGH25	AGHS4	-0.2590	0.2206	161.5	-1.17	0.2420	0.8488
Year First	AGHS1	AGHS2	0.01642	0.2206	161.5	0.07	0.9408	1.0000
Year First	AGHS1	AGHS3	0.1042	0.2206	161.5	0.47	0.6374	0.9971
Year First	AGHS1	AGHS4	0.02282	0.2206	161.5	0.10	0.9177	1.0000
Year First	AGHS2	AGHS3	0.08776	0.2206	161.5	0.40	0.6913	0.9987
Year First	AGHS2	AGHS4	0.006401	0.2206	161.5	0.03	0.9769	1.0000
Year First	AGHS3	AGHS4	-0.08136	0.2206	161.5	-0.37	0.7127	0.9991

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Family	Estimate
Year First	AGHS1	0.9817	A
Year First			A
Year First	AGHS2	0.9653	A
Year First			A
Year First	AGHS4	0.9589	A
Year First			A
Year First	AGHS3	0.8775	A
Year First			A
Year First	AGH25	0.6998	A
Year First			A
Year First	AGH2	0.6970	A

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Fourth	5	161.5	20.56	<.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Family	_Family	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Fourth	AGH2	AGH25	-0.3428	0.2206	161.5	-1.55	0.1221	0.6294
Year Fourth	AGH2	AGHS1	-1.5691	0.2206	161.5	-7.11	<.0001	<.0001
Year Fourth	AGH2	AGHS2	-1.7651	0.2206	161.5	-8.00	<.0001	<.0001
Year Fourth	AGH2	AGHS3	-0.9453	0.2206	161.5	-4.29	<.0001	0.0004
Year Fourth	AGH2	AGHS4	-1.3959	0.2206	161.5	-6.33	<.0001	<.0001
Year Fourth	AGH25	AGHS1	-1.2262	0.2206	161.5	-5.56	<.0001	<.0001
Year Fourth	AGH25	AGHS2	-1.4223	0.2206	161.5	-6.45	<.0001	<.0001
Year Fourth	AGH25	AGHS3	-0.6025	0.2206	161.5	-2.73	0.0070	0.0727
Year Fourth	AGH25	AGHS4	-1.0531	0.2206	161.5	-4.77	<.0001	<.0001
Year Fourth	AGHS1	AGHS2	-0.1960	0.2206	161.5	-0.89	0.3754	0.9490
Year Fourth	AGHS1	AGHS3	0.6237	0.2206	161.5	2.83	0.0053	0.0564
Year Fourth	AGHS1	AGHS4	0.1731	0.2206	161.5	0.78	0.4337	0.9699
Year Fourth	AGHS2	AGHS3	0.8198	0.2206	161.5	3.72	0.0003	0.0034
Year Fourth	AGHS2	AGHS4	0.3692	0.2206	161.5	1.67	0.0961	0.5505
Year Fourth	AGHS3	AGHS4	-0.4506	0.2206	161.5	-2.04	0.0427	0.3212

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Family	Estimate
Year Fourth	AGHS2	5.5407		A
Year Fourth				A
Year Fourth	AGHS1	5.3447	B	A
Year Fourth			B	A
Year Fourth	AGHS4	5.1715	B	A
Year Fourth			B
Year Fourth	AGHS3	4.7209	B	C
Year Fourth				C
Year Fourth	AGH25	4.1184	D	C
Year Fourth			D
Year Fourth	AGH2	3.7756	D

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Initial	5	161.5	0.14	0.9815

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Family	_Family	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Initial	AGH2	AGH25	0.02700	0.2206	161.5	0.12	0.9027	1.0000
Year Initial	AGH2	AGHS1	-0.09858	0.2206	161.5	-0.45	0.6555	0.9977
Year Initial	AGH2	AGHS2	-0.02001	0.2206	161.5	-0.09	0.9278	1.0000
Year Initial	AGH2	AGHS3	-0.09166	0.2206	161.5	-0.42	0.6783	0.9984
Year Initial	AGH2	AGHS4	-0.1142	0.2206	161.5	-0.52	0.6053	0.9955
Year Initial	AGH25	AGHS1	-0.1256	0.2206	161.5	-0.57	0.5699	0.9929
Year Initial	AGH25	AGHS2	-0.04701	0.2206	161.5	-0.21	0.8315	0.9999
Year Initial	AGH25	AGHS3	-0.1187	0.2206	161.5	-0.54	0.5914	0.9946
Year Initial	AGH25	AGHS4	-0.1412	0.2206	161.5	-0.64	0.5230	0.9878
Year Initial	AGHS1	AGHS2	0.07857	0.2206	161.5	0.36	0.7221	0.9992
Year Initial	AGHS1	AGHS3	0.006927	0.2206	161.5	0.03	0.9750	1.0000
Year Initial	AGHS1	AGHS4	-0.01562	0.2206	161.5	-0.07	0.9436	1.0000
Year Initial	AGHS2	AGHS3	-0.07165	0.2206	161.5	-0.32	0.7457	0.9995
Year Initial	AGHS2	AGHS4	-0.09420	0.2206	161.5	-0.43	0.6699	0.9982
Year Initial	AGHS3	AGHS4	-0.02255	0.2206	161.5	-0.10	0.9187	1.0000

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Family	Estimate
Year Initial	AGHS4	0.5520	A
Year Initial			A
Year Initial	AGHS1	0.5363	A
Year Initial			A
Year Initial	AGHS3	0.5294	A
Year Initial			A
Year Initial	AGHS2	0.4578	A
Year Initial			A
Year Initial	AGH2	0.4378	A
Year Initial			A
Year Initial	AGH25	0.4108	A

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Second	5	161.5	5.65	<.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Family	_Family	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Second	AGH2	AGH25	-0.09287	0.2206	161.5	-0.42	0.6743	0.9983
Year Second	AGH2	AGHS1	-0.7712	0.2206	161.5	-3.50	0.0006	0.0073
Year Second	AGH2	AGHS2	-0.9229	0.2206	161.5	-4.18	<.0001	0.0006
Year Second	AGH2	AGHS3	-0.5018	0.2206	161.5	-2.28	0.0242	0.2082
Year Second	AGH2	AGHS4	-0.6426	0.2206	161.5	-2.91	0.0041	0.0446
Year Second	AGH25	AGHS1	-0.6784	0.2206	161.5	-3.08	0.0025	0.0279
Year Second	AGH25	AGHS2	-0.8300	0.2206	161.5	-3.76	0.0002	0.0028
Year Second	AGH25	AGHS3	-0.4090	0.2206	161.5	-1.85	0.0656	0.4332
Year Second	AGH25	AGHS4	-0.5497	0.2206	161.5	-2.49	0.0137	0.1303
Year Second	AGHS1	AGHS2	-0.1517	0.2206	161.5	-0.69	0.4927	0.9832
Year Second	AGHS1	AGHS3	0.2694	0.2206	161.5	1.22	0.2237	0.8261
Year Second	AGHS1	AGHS4	0.1286	0.2206	161.5	0.58	0.5606	0.9921
Year Second	AGHS2	AGHS3	0.4211	0.2206	161.5	1.91	0.0580	0.3991
Year Second	AGHS2	AGHS4	0.2803	0.2206	161.5	1.27	0.2056	0.8007
Year Second	AGHS3	AGHS4	-0.1408	0.2206	161.5	-0.64	0.5243	0.9880

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Family	Estimate
Year Second	AGHS2	2.5035		A
Year Second				A
Year Second	AGHS1	2.3519		A
Year Second				A
Year Second	AGHS4	2.2232	B	A
Year Second			B	A
Year Second	AGHS3	2.0825	B	A	C
Year Second			B		C
Year Second	AGH25	1.6735	B		C
Year Second					C
Year Second	AGH2	1.5806			C

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Seventh	5	161.5	15.75	<.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Family	_Family	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Seventh	AGH2	AGH25	0.1012	0.2206	161.5	0.46	0.6469	0.9974
Year Seventh	AGH2	AGHS1	-1.2296	0.2206	161.5	-5.57	<.0001	<.0001
Year Seventh	AGH2	AGHS2	-1.2242	0.2206	161.5	-5.55	<.0001	<.0001
Year Seventh	AGH2	AGHS3	-0.5381	0.2206	161.5	-2.44	0.0158	0.1468
Year Seventh	AGH2	AGHS4	-1.1547	0.2206	161.5	-5.23	<.0001	<.0001
Year Seventh	AGH25	AGHS1	-1.3308	0.2206	161.5	-6.03	<.0001	<.0001
Year Seventh	AGH25	AGHS2	-1.3254	0.2206	161.5	-6.01	<.0001	<.0001
Year Seventh	AGH25	AGHS3	-0.6393	0.2206	161.5	-2.90	0.0043	0.0465
Year Seventh	AGH25	AGHS4	-1.2559	0.2206	161.5	-5.69	<.0001	<.0001
Year Seventh	AGHS1	AGHS2	0.005355	0.2206	161.5	0.02	0.9807	1.0000
Year Seventh	AGHS1	AGHS3	0.6915	0.2206	161.5	3.13	0.0020	0.0234
Year Seventh	AGHS1	AGHS4	0.07490	0.2206	161.5	0.34	0.7346	0.9994
Year Seventh	AGHS2	AGHS3	0.6861	0.2206	161.5	3.11	0.0022	0.0251
Year Seventh	AGHS2	AGHS4	0.06954	0.2206	161.5	0.32	0.7530	0.9996
Year Seventh	AGHS3	AGHS4	-0.6166	0.2206	161.5	-2.80	0.0058	0.0615

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Family	Estimate
Year Seventh	AGHS1	9.6536		A
Year Seventh				A
Year Seventh	AGHS2	9.6482		A
Year Seventh				A
Year Seventh	AGHS4	9.5787	B	A
Year Seventh			B
Year Seventh	AGHS3	8.9621	B	C
Year Seventh				C
Year Seventh	AGH2	8.4240	D	C
Year Seventh			D
Year Seventh	AGH25	8.3228	D

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Sixth	5	161.5	27.96	<.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Family	_Family	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Sixth	AGH2	AGH25	-0.3353	0.2206	161.5	-1.52	0.1304	0.6515
Year Sixth	AGH2	AGHS1	-1.8840	0.2206	161.5	-8.54	<.0001	<.0001
Year Sixth	AGH2	AGHS2	-1.9521	0.2206	161.5	-8.85	<.0001	<.0001
Year Sixth	AGH2	AGHS3	-1.2217	0.2206	161.5	-5.54	<.0001	<.0001
Year Sixth	AGH2	AGHS4	-1.6303	0.2206	161.5	-7.39	<.0001	<.0001
Year Sixth	AGH25	AGHS1	-1.5486	0.2206	161.5	-7.02	<.0001	<.0001
Year Sixth	AGH25	AGHS2	-1.6167	0.2206	161.5	-7.33	<.0001	<.0001
Year Sixth	AGH25	AGHS3	-0.8864	0.2206	161.5	-4.02	<.0001	0.0011
Year Sixth	AGH25	AGHS4	-1.2950	0.2206	161.5	-5.87	<.0001	<.0001
Year Sixth	AGHS1	AGHS2	-0.06811	0.2206	161.5	-0.31	0.7579	0.9996
Year Sixth	AGHS1	AGHS3	0.6623	0.2206	161.5	3.00	0.0031	0.0346
Year Sixth	AGHS1	AGHS4	0.2536	0.2206	161.5	1.15	0.2519	0.8599
Year Sixth	AGHS2	AGHS3	0.7304	0.2206	161.5	3.31	0.0011	0.0135
Year Sixth	AGHS2	AGHS4	0.3217	0.2206	161.5	1.46	0.1466	0.6909
Year Sixth	AGHS3	AGHS4	-0.4086	0.2206	161.5	-1.85	0.0658	0.4341

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Family	Estimate
Year Sixth	AGHS2	8.1945		A
Year Sixth				A
Year Sixth	AGHS1	8.1264		A
Year Sixth				A
Year Sixth	AGHS4	7.8728	B	A
Year Sixth			B
Year Sixth	AGHS3	7.4642	B
Year Sixth
Year Sixth	AGH25	6.5778		C
Year Sixth				C
Year Sixth	AGH2	6.2425		C

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Third	5	161.5	8.12	<.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Family	_Family	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Third	AGH2	AGH25	-0.1341	0.2206	161.5	-0.61	0.5442	0.9904
Year Third	AGH2	AGHS1	-1.0244	0.2206	161.5	-4.64	<.0001	<.0001
Year Third	AGH2	AGHS2	-1.0013	0.2206	161.5	-4.54	<.0001	0.0001
Year Third	AGH2	AGHS3	-0.3330	0.2206	161.5	-1.51	0.1331	0.6585
Year Third	AGH2	AGHS4	-0.7438	0.2206	161.5	-3.37	0.0009	0.0110
Year Third	AGH25	AGHS1	-0.8903	0.2206	161.5	-4.04	<.0001	0.0010
Year Third	AGH25	AGHS2	-0.8673	0.2206	161.5	-3.93	0.0001	0.0015
Year Third	AGH25	AGHS3	-0.1989	0.2206	161.5	-0.90	0.3686	0.9458
Year Third	AGH25	AGHS4	-0.6097	0.2206	161.5	-2.76	0.0064	0.0667
Year Third	AGHS1	AGHS2	0.02302	0.2206	161.5	0.10	0.9170	1.0000
Year Third	AGHS1	AGHS3	0.6914	0.2206	161.5	3.13	0.0020	0.0234
Year Third	AGHS1	AGHS4	0.2806	0.2206	161.5	1.27	0.2052	0.8000
Year Third	AGHS2	AGHS3	0.6684	0.2206	161.5	3.03	0.0028	0.0319
Year Third	AGHS2	AGHS4	0.2576	0.2206	161.5	1.17	0.2447	0.8518
Year Third	AGHS3	AGHS4	-0.4108	0.2206	161.5	-1.86	0.0643	0.4278

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Family	Estimate
Year Third	AGHS1	3.5000		A
Year Third				A
Year Third	AGHS2	3.4770		A
Year Third				A
Year Third	AGHS4	3.2194	B	A
Year Third			B
Year Third	AGHS3	2.8086	B	C
Year Third			B	C
Year Third	AGH25	2.6097	B	C
Year Third				C
Year Third	AGH2	2.4756		C

LS-Means Results

Family LS-Means

Family Least Squares Means
Family	Estimate	Standard Error	DF	t Value	Pr > |t|
AGH2	3.6292	0.1008	36	36.01	<.0001
AGH25	3.7903	0.1008	36	37.60	<.0001
AGHS1	4.6630	0.1008	36	46.26	<.0001
AGHS2	4.7374	0.1008	36	47.00	<.0001
AGHS3	4.1983	0.1008	36	41.65	<.0001
AGHS4	4.5520	0.1008	36	45.16	<.0001

Family Diffs

Differences of Family Least Squares Means Adjustment for Multiple Comparisons: Tukey
Family	_Family	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
AGH2	AGH25	-0.1611	0.1425	36	-1.13	0.2660	0.8656
AGH2	AGHS1	-1.0338	0.1425	36	-7.25	<.0001	<.0001
AGH2	AGHS2	-1.1081	0.1425	36	-7.77	<.0001	<.0001
AGH2	AGHS3	-0.5691	0.1425	36	-3.99	0.0003	0.0039
AGH2	AGHS4	-0.9227	0.1425	36	-6.47	<.0001	<.0001
AGH25	AGHS1	-0.8727	0.1425	36	-6.12	<.0001	<.0001
AGH25	AGHS2	-0.9471	0.1425	36	-6.64	<.0001	<.0001
AGH25	AGHS3	-0.4080	0.1425	36	-2.86	0.0070	0.0700
AGH25	AGHS4	-0.7616	0.1425	36	-5.34	<.0001	<.0001
AGHS1	AGHS2	-0.07436	0.1425	36	-0.52	0.6051	0.9950
AGHS1	AGHS3	0.4647	0.1425	36	3.26	0.0024	0.0272
AGHS1	AGHS4	0.1111	0.1425	36	0.78	0.4410	0.9693
AGHS2	AGHS3	0.5390	0.1425	36	3.78	0.0006	0.0069
AGHS2	AGHS4	0.1854	0.1425	36	1.30	0.2016	0.7827
AGHS3	AGHS4	-0.3536	0.1425	36	-2.48	0.0179	0.1569

Family Lines

Tukey Grouping for Family Least Squares Means (Alpha=0.05)
LS-means with the same letter are not significantly different.
Family	Estimate
AGHS2	4.7374		A
			A
AGHS1	4.6630		A
			A
AGHS4	4.5520	B	A
		B
AGHS3	4.1983	B	C
			C
AGH25	3.7903	D	C
		D
AGH2	3.6292	D

SAS Output Crown

htmltools::includeHTML("https://raw.githubusercontent.com/azd169/SAS_outputs/main/hane-results-crown.html")

## Warning: `includeHTML()` was provided a `path` that appears to be a complete HTML document.
## ✖ Path: https://raw.githubusercontent.com/azd169/SAS_outputs/main/hane-results-crown.html
## ℹ Use `tags$iframe()` to include an HTML document. You can either ensure `path` is accessible in your app or document (see e.g. `shiny::addResourcePath()`) and pass the relative path to the `src` argument. Or you can read the contents of `path` and pass the contents to `srcdoc`.

Results: hane.sas

The GLIMMIX Procedure

The Glimmix Procedure

Model Information

Model Information
Data Set	WORK.CROWN
Response Variable	Crown
Response Distribution	Gaussian
Link Function	Identity
Variance Function	Default
Variance Matrix Blocked By	GLB
Estimation Technique	Restricted Maximum Likelihood
Degrees of Freedom Method	Kenward-Roger
Fixed Effects SE Adjustment	Kenward-Roger

Class Level Information

Class Level Information
Class	Levels	Values
Year	3	Fifth Seventh Sixth
Location	2	Hill Farm LA Tech
Block	4	1 2 3 4
Family	6	AGH2 AGH25 AGHS1 AGHS2 AGHS3 AGHS4
GLB	48	AGH25_Hill Farm_1 AGH25_Hill Farm_2 AGH25_Hill Farm_3 AGH25_Hill Farm_4 AGH25_LA Tech_1 AGH25_LA Tech_2 AGH25_LA Tech_3 AGH25_LA Tech_4 AGH2_Hill Farm_1 AGH2_Hill Farm_2 AGH2_Hill Farm_3 AGH2_Hill Farm_4 AGH2_LA Tech_1 AGH2_LA Tech_2 AGH2_LA Tech_3 AGH2_LA Tech_4 AGHS1_Hill Farm_1 AGHS1_Hill Farm_2 AGHS1_Hill Farm_3 AGHS1_Hill Farm_4 AGHS1_LA Tech_1 AGHS1_LA Tech_2 AGHS1_LA Tech_3 AGHS1_LA Tech_4 AGHS2_Hill Farm_1 AGHS2_Hill Farm_2 AGHS2_Hill Farm_3 AGHS2_Hill Farm_4 AGHS2_LA Tech_1 AGHS2_LA Tech_2 AGHS2_LA Tech_3 AGHS2_LA Tech_4 AGHS3_Hill Farm_1 AGHS3_Hill Farm_2 AGHS3_Hill Farm_3 AGHS3_Hill Farm_4 AGHS3_LA Tech_1 AGHS3_LA Tech_2 AGHS3_LA Tech_3 AGHS3_LA Tech_4 AGHS4_Hill Farm_1 AGHS4_Hill Farm_2 AGHS4_Hill Farm_3 AGHS4_Hill Farm_4 AGHS4_LA Tech_1 AGHS4_LA Tech_2 AGHS4_LA Tech_3 AGHS4_LA Tech_4

Number of Observations

Number of Observations Read	144
Number of Observations Used	144

Dimensions

Dimensions
G-side Cov. Parameters	2
R-side Cov. Parameters	1
Columns in X	84
Columns in Z per Subject	1
Subjects (Blocks in V)	48
Max Obs per Subject	3

Optimization Information

Optimization Information
Optimization Technique	Dual Quasi-Newton
Parameters in Optimization	2
Lower Boundaries	2
Upper Boundaries	1
Fixed Effects	Profiled
Residual Variance	Profiled
Starting From	Data

Iteration History

Iteration History
Iteration	Restarts	Evaluations	Objective Function	Change	Max Gradient
0	0	4	327.11392315	.	7.57E-13

Convergence Status

Convergence criterion (ABSGCONV=0.00001) satisfied.

Hessian Matrix

Hessian Matrix
Cov Parm	Subject	Row	Col1	Col2
Variance	GLB	1	14.9172
AR(1)	GLB	2

Fit Statistics

Fit Statistics
-2 Res Log Likelihood	327.11
AIC (smaller is better)	333.11
AICC (smaller is better)	333.34
BIC (smaller is better)	338.73
CAIC (smaller is better)	341.73
HQIC (smaller is better)	335.24
Generalized Chi-Square	52.75
Gener. Chi-Square / DF	0.49

Covariance Parameter Estimates

Covariance Parameter Estimates
Cov Parm	Subject	Estimate	Standard Error
Variance	GLB	0.4567	0.1485
AR(1)	GLB	0	.
Residual		0.4884	0.08140

Solutions for Fixed Effects

Solutions for Fixed Effects
Effect	Year	Location	Family	Estimate	Standard Error	DF	t Value	Pr > |t|
Intercept				4.6956	0.4861	73.62	9.66	<.0001
Year	Fifth			-3.0771	0.4942	72	-6.23	<.0001
Year	Seventh			1.0079	0.4942	72	2.04	0.0451
Year	Sixth			0	.	.	.	.
Location		Hill Farm		-2.5922	0.6874	73.62	-3.77	0.0003
Location		LA Tech		0	.	.	.	.
Year*Location	Fifth	Hill Farm		1.4553	0.6989	72	2.08	0.0409
Year*Location	Fifth	LA Tech		0	.	.	.	.
Year*Location	Seventh	Hill Farm		-0.1176	0.6989	72	-0.17	0.8668
Year*Location	Seventh	LA Tech		0	.	.	.	.
Year*Location	Sixth	Hill Farm		0	.	.	.	.
Year*Location	Sixth	LA Tech		0	.	.	.	.
Family			AGH2	-3.0492	0.6874	73.62	-4.44	<.0001
Family			AGH25	-2.3896	0.6874	73.62	-3.48	0.0009
Family			AGHS1	0.5941	0.6874	73.62	0.86	0.3902
Family			AGHS2	0.8519	0.6874	73.62	1.24	0.2192
Family			AGHS3	-0.08741	0.6874	73.62	-0.13	0.8992
Family			AGHS4	0	.	.	.	.
Year*Family	Fifth		AGH2	2.6497	0.6989	72	3.79	0.0003
Year*Family	Fifth		AGH25	1.7974	0.6989	72	2.57	0.0122
Year*Family	Fifth		AGHS1	-0.1967	0.6989	72	-0.28	0.7791
Year*Family	Fifth		AGHS2	0.2015	0.6989	72	0.29	0.7739
Year*Family	Fifth		AGHS3	0.3601	0.6989	72	0.52	0.6079
Year*Family	Fifth		AGHS4	0	.	.	.	.
Year*Family	Seventh		AGH2	0.2702	0.6989	72	0.39	0.7002
Year*Family	Seventh		AGH25	-0.4015	0.6989	72	-0.57	0.5674
Year*Family	Seventh		AGHS1	-0.2398	0.6989	72	-0.34	0.7325
Year*Family	Seventh		AGHS2	-0.3092	0.6989	72	-0.44	0.6595
Year*Family	Seventh		AGHS3	-0.6890	0.6989	72	-0.99	0.3275
Year*Family	Seventh		AGHS4	0	.	.	.	.
Year*Family	Sixth		AGH2	0	.	.	.	.
Year*Family	Sixth		AGH25	0	.	.	.	.
Year*Family	Sixth		AGHS1	0	.	.	.	.
Year*Family	Sixth		AGHS2	0	.	.	.	.
Year*Family	Sixth		AGHS3	0	.	.	.	.
Year*Family	Sixth		AGHS4	0	.	.	.	.
Location*Family		Hill Farm	AGH2	1.8355	0.9722	73.62	1.89	0.0630
Location*Family		Hill Farm	AGH25	1.4387	0.9722	73.62	1.48	0.1432
Location*Family		Hill Farm	AGHS1	0.1333	0.9722	73.62	0.14	0.8913
Location*Family		Hill Farm	AGHS2	0.4140	0.9722	73.62	0.43	0.6714
Location*Family		Hill Farm	AGHS3	0.1992	0.9722	73.62	0.20	0.8382
Location*Family		Hill Farm	AGHS4	0	.	.	.	.
Location*Family		LA Tech	AGH2	0	.	.	.	.
Location*Family		LA Tech	AGH25	0	.	.	.	.
Location*Family		LA Tech	AGHS1	0	.	.	.	.
Location*Family		LA Tech	AGHS2	0	.	.	.	.
Location*Family		LA Tech	AGHS3	0	.	.	.	.
Location*Family		LA Tech	AGHS4	0	.	.	.	.
Year*Location*Family	Fifth	Hill Farm	AGH2	-1.4704	0.9883	72	-1.49	0.1412
Year*Location*Family	Fifth	Hill Farm	AGH25	-0.8170	0.9883	72	-0.83	0.4112
Year*Location*Family	Fifth	Hill Farm	AGHS1	-0.1838	0.9883	72	-0.19	0.8530
Year*Location*Family	Fifth	Hill Farm	AGHS2	-1.0488	0.9883	72	-1.06	0.2921
Year*Location*Family	Fifth	Hill Farm	AGHS3	-0.4048	0.9883	72	-0.41	0.6833
Year*Location*Family	Fifth	Hill Farm	AGHS4	0	.	.	.	.
Year*Location*Family	Fifth	LA Tech	AGH2	0	.	.	.	.
Year*Location*Family	Fifth	LA Tech	AGH25	0	.	.	.	.
Year*Location*Family	Fifth	LA Tech	AGHS1	0	.	.	.	.
Year*Location*Family	Fifth	LA Tech	AGHS2	0	.	.	.	.
Year*Location*Family	Fifth	LA Tech	AGHS3	0	.	.	.	.
Year*Location*Family	Fifth	LA Tech	AGHS4	0	.	.	.	.
Year*Location*Family	Seventh	Hill Farm	AGH2	-0.5397	0.9883	72	-0.55	0.5867
Year*Location*Family	Seventh	Hill Farm	AGH25	-0.4144	0.9883	72	-0.42	0.6762
Year*Location*Family	Seventh	Hill Farm	AGHS1	0.1666	0.9883	72	0.17	0.8666
Year*Location*Family	Seventh	Hill Farm	AGHS2	0.3237	0.9883	72	0.33	0.7443
Year*Location*Family	Seventh	Hill Farm	AGHS3	0.5593	0.9883	72	0.57	0.5732
Year*Location*Family	Seventh	Hill Farm	AGHS4	0	.	.	.	.
Year*Location*Family	Seventh	LA Tech	AGH2	0	.	.	.	.
Year*Location*Family	Seventh	LA Tech	AGH25	0	.	.	.	.
Year*Location*Family	Seventh	LA Tech	AGHS1	0	.	.	.	.
Year*Location*Family	Seventh	LA Tech	AGHS2	0	.	.	.	.
Year*Location*Family	Seventh	LA Tech	AGHS3	0	.	.	.	.
Year*Location*Family	Seventh	LA Tech	AGHS4	0	.	.	.	.
Year*Location*Family	Sixth	Hill Farm	AGH2	0	.	.	.	.
Year*Location*Family	Sixth	Hill Farm	AGH25	0	.	.	.	.
Year*Location*Family	Sixth	Hill Farm	AGHS1	0	.	.	.	.
Year*Location*Family	Sixth	Hill Farm	AGHS2	0	.	.	.	.
Year*Location*Family	Sixth	Hill Farm	AGHS3	0	.	.	.	.
Year*Location*Family	Sixth	Hill Farm	AGHS4	0	.	.	.	.
Year*Location*Family	Sixth	LA Tech	AGH2	0	.	.	.	.
Year*Location*Family	Sixth	LA Tech	AGH25	0	.	.	.	.
Year*Location*Family	Sixth	LA Tech	AGHS1	0	.	.	.	.
Year*Location*Family	Sixth	LA Tech	AGHS2	0	.	.	.	.
Year*Location*Family	Sixth	LA Tech	AGHS3	0	.	.	.	.
Year*Location*Family	Sixth	LA Tech	AGHS4	0	.	.	.	.

Type III Tests of Fixed Effects

Type III Tests of Fixed Effects
Effect	Num DF	Den DF	F Value	Pr > F
Year	2	72	177.27	<.0001
Location	1	36	55.25	<.0001
Year*Location	2	72	6.01	0.0039
Family	5	36	12.64	<.0001
Year*Family	10	72	5.39	<.0001
Location*Family	5	36	0.82	0.5439
Year*Location*Family	10	72	0.51	0.8776

Solution for Random Effects

Solution for Random Effects
Effect	Subject	Estimate	Std Err Pred	DF	t Value	Pr > |t|	Alpha	Lower	Upper
Intercept	GLB AGH25_Hill Farm_1	-0.1422	0.4578	42.76	-0.31	0.7575	0.05	-1.0655	0.7811
Intercept	GLB AGH25_Hill Farm_2	0.2869	0.4578	42.76	0.63	0.5341	0.05	-0.6364	1.2102
Intercept	GLB AGH25_Hill Farm_3	-0.2848	0.4578	42.76	-0.62	0.5371	0.05	-1.2081	0.6385
Intercept	GLB AGH25_Hill Farm_4	0.1401	0.4578	42.76	0.31	0.7611	0.05	-0.7832	1.0634
Intercept	GLB AGH25_LA Tech_1	-0.6896	0.4578	42.76	-1.51	0.1393	0.05	-1.6129	0.2337
Intercept	GLB AGH25_LA Tech_2	0.1718	0.4578	42.76	0.38	0.7093	0.05	-0.7515	1.0951
Intercept	GLB AGH25_LA Tech_3	0.2293	0.4578	42.76	0.50	0.6190	0.05	-0.6940	1.1526
Intercept	GLB AGH25_LA Tech_4	0.2885	0.4578	42.76	0.63	0.5319	0.05	-0.6348	1.2118
Intercept	GLB AGH2_Hill Farm_1	-0.3502	0.4578	42.76	-0.77	0.4484	0.05	-1.2735	0.5731
Intercept	GLB AGH2_Hill Farm_2	0.2129	0.4578	42.76	0.47	0.6442	0.05	-0.7104	1.1362
Intercept	GLB AGH2_Hill Farm_3	0.2285	0.4578	42.76	0.50	0.6201	0.05	-0.6948	1.1518
Intercept	GLB AGH2_Hill Farm_4	-0.09127	0.4578	42.76	-0.20	0.8429	0.05	-1.0146	0.8320
Intercept	GLB AGH2_LA Tech_1	-0.5441	0.4578	42.76	-1.19	0.2412	0.05	-1.4673	0.3792
Intercept	GLB AGH2_LA Tech_2	-0.2853	0.4578	42.76	-0.62	0.5364	0.05	-1.2086	0.6380
Intercept	GLB AGH2_LA Tech_3	0.8694	0.4578	42.76	1.90	0.0643	0.05	-0.05387	1.7927
Intercept	GLB AGH2_LA Tech_4	-0.04006	0.4578	42.76	-0.09	0.9307	0.05	-0.9634	0.8832
Intercept	GLB AGHS1_Hill Farm_1	-0.4909	0.4578	42.76	-1.07	0.2896	0.05	-1.4142	0.4324
Intercept	GLB AGHS1_Hill Farm_2	0.5528	0.4578	42.76	1.21	0.2338	0.05	-0.3705	1.4761
Intercept	GLB AGHS1_Hill Farm_3	0.4166	0.4578	42.76	0.91	0.3679	0.05	-0.5067	1.3399
Intercept	GLB AGHS1_Hill Farm_4	-0.4785	0.4578	42.76	-1.05	0.3017	0.05	-1.4018	0.4448
Intercept	GLB AGHS1_LA Tech_1	-1.1449	0.4578	42.76	-2.50	0.0163	0.05	-2.0682	-0.2216
Intercept	GLB AGHS1_LA Tech_2	-0.2571	0.4578	42.76	-0.56	0.5772	0.05	-1.1804	0.6662
Intercept	GLB AGHS1_LA Tech_3	0.8501	0.4578	42.76	1.86	0.0702	0.05	-0.07321	1.7734
Intercept	GLB AGHS1_LA Tech_4	0.5519	0.4578	42.76	1.21	0.2346	0.05	-0.3714	1.4752
Intercept	GLB AGHS2_Hill Farm_1	-0.9490	0.4578	42.76	-2.07	0.0442	0.05	-1.8723	-0.02572
Intercept	GLB AGHS2_Hill Farm_2	0.3246	0.4578	42.76	0.71	0.4821	0.05	-0.5987	1.2479
Intercept	GLB AGHS2_Hill Farm_3	0.2169	0.4578	42.76	0.47	0.6381	0.05	-0.7064	1.1402
Intercept	GLB AGHS2_Hill Farm_4	0.4075	0.4578	42.76	0.89	0.3783	0.05	-0.5158	1.3308
Intercept	GLB AGHS2_LA Tech_1	-0.9095	0.4578	42.76	-1.99	0.0534	0.05	-1.8328	0.01377
Intercept	GLB AGHS2_LA Tech_2	0.06896	0.4578	42.76	0.15	0.8810	0.05	-0.8543	0.9923
Intercept	GLB AGHS2_LA Tech_3	0.7353	0.4578	42.76	1.61	0.1156	0.05	-0.1880	1.6586
Intercept	GLB AGHS2_LA Tech_4	0.1053	0.4578	42.76	0.23	0.8192	0.05	-0.8180	1.0286
Intercept	GLB AGHS3_Hill Farm_1	-0.4208	0.4578	42.76	-0.92	0.3631	0.05	-1.3441	0.5025
Intercept	GLB AGHS3_Hill Farm_2	-0.4823	0.4578	42.76	-1.05	0.2980	0.05	-1.4056	0.4410
Intercept	GLB AGHS3_Hill Farm_3	0.6032	0.4578	42.76	1.32	0.1946	0.05	-0.3201	1.5265
Intercept	GLB AGHS3_Hill Farm_4	0.2999	0.4578	42.76	0.66	0.5159	0.05	-0.6234	1.2232
Intercept	GLB AGHS3_LA Tech_1	-0.6954	0.4578	42.76	-1.52	0.1361	0.05	-1.6187	0.2279
Intercept	GLB AGHS3_LA Tech_2	-0.3490	0.4578	42.76	-0.76	0.4500	0.05	-1.2723	0.5743
Intercept	GLB AGHS3_LA Tech_3	0.9425	0.4578	42.76	2.06	0.0456	0.05	0.01919	1.8658
Intercept	GLB AGHS3_LA Tech_4	0.1019	0.4578	42.76	0.22	0.8249	0.05	-0.8214	1.0252
Intercept	GLB AGHS4_Hill Farm_1	-0.7912	0.4578	42.76	-1.73	0.0911	0.05	-1.7145	0.1321
Intercept	GLB AGHS4_Hill Farm_2	0.06875	0.4578	42.76	0.15	0.8813	0.05	-0.8545	0.9920
Intercept	GLB AGHS4_Hill Farm_3	0.005591	0.4578	42.76	0.01	0.9903	0.05	-0.9177	0.9289
Intercept	GLB AGHS4_Hill Farm_4	0.7168	0.4578	42.76	1.57	0.1247	0.05	-0.2065	1.6401
Intercept	GLB AGHS4_LA Tech_1	-0.5740	0.4578	42.76	-1.25	0.2167	0.05	-1.4973	0.3493
Intercept	GLB AGHS4_LA Tech_2	0.2233	0.4578	42.76	0.49	0.6281	0.05	-0.7000	1.1466
Intercept	GLB AGHS4_LA Tech_3	0.2200	0.4578	42.76	0.48	0.6332	0.05	-0.7033	1.1433
Intercept	GLB AGHS4_LA Tech_4	0.1306	0.4578	42.76	0.29	0.7768	0.05	-0.7927	1.0539

Slices

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Fifth	5	73.62	1.24	0.2982

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Family	_Family	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Fifth	AGH2	AGH25	0.06438	0.4861	73.62	0.13	0.8950	1.0000
Year Fifth	AGH2	AGHS1	-0.5890	0.4861	73.62	-1.21	0.2295	0.8297
Year Fifth	AGH2	AGHS2	-0.9529	0.4861	73.62	-1.96	0.0537	0.3750
Year Fifth	AGH2	AGHS3	-0.3868	0.4861	73.62	-0.80	0.4287	0.9674
Year Fifth	AGH2	AGHS4	-0.2169	0.4861	73.62	-0.45	0.6567	0.9977
Year Fifth	AGH25	AGHS1	-0.6534	0.4861	73.62	-1.34	0.1830	0.7594
Year Fifth	AGH25	AGHS2	-1.0173	0.4861	73.62	-2.09	0.0398	0.3027
Year Fifth	AGH25	AGHS3	-0.4512	0.4861	73.62	-0.93	0.3563	0.9379
Year Fifth	AGH25	AGHS4	-0.2813	0.4861	73.62	-0.58	0.5646	0.9921
Year Fifth	AGHS1	AGHS2	-0.3639	0.4861	73.62	-0.75	0.4565	0.9750
Year Fifth	AGHS1	AGHS3	0.2023	0.4861	73.62	0.42	0.6786	0.9983
Year Fifth	AGHS1	AGHS4	0.3721	0.4861	73.62	0.77	0.4464	0.9724
Year Fifth	AGHS2	AGHS3	0.5661	0.4861	73.62	1.16	0.2479	0.8520
Year Fifth	AGHS2	AGHS4	0.7360	0.4861	73.62	1.51	0.1343	0.6564
Year Fifth	AGHS3	AGHS4	0.1699	0.4861	73.62	0.35	0.7277	0.9993

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Family	Estimate
Year Fifth	AGHS2	1.7861	A
Year Fifth			A
Year Fifth	AGHS1	1.4222	A
Year Fifth			A
Year Fifth	AGHS3	1.2199	A
Year Fifth			A
Year Fifth	AGHS4	1.0501	A
Year Fifth			A
Year Fifth	AGH2	0.8332	A
Year Fifth			A
Year Fifth	AGH25	0.7688	A

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Seventh	5	73.62	15.34	<.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Family	_Family	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Seventh	AGH2	AGH25	0.1479	0.4861	73.62	0.30	0.7619	0.9996
Year Seventh	AGH2	AGHS1	-2.6354	0.4861	73.62	-5.42	<.0001	<.0001
Year Seventh	AGH2	AGHS2	-3.0426	0.4861	73.62	-6.26	<.0001	<.0001
Year Seventh	AGH2	AGHS3	-1.7339	0.4861	73.62	-3.57	0.0006	0.0082
Year Seventh	AGH2	AGHS4	-2.1311	0.4861	73.62	-4.38	<.0001	0.0005
Year Seventh	AGH25	AGHS1	-2.7833	0.4861	73.62	-5.73	<.0001	<.0001
Year Seventh	AGH25	AGHS2	-3.1904	0.4861	73.62	-6.56	<.0001	<.0001
Year Seventh	AGH25	AGHS3	-1.8818	0.4861	73.62	-3.87	0.0002	0.0031
Year Seventh	AGH25	AGHS4	-2.2789	0.4861	73.62	-4.69	<.0001	0.0002
Year Seventh	AGHS1	AGHS2	-0.4072	0.4861	73.62	-0.84	0.4049	0.9594
Year Seventh	AGHS1	AGHS3	0.9015	0.4861	73.62	1.85	0.0677	0.4382
Year Seventh	AGHS1	AGHS4	0.5043	0.4861	73.62	1.04	0.3029	0.9036
Year Seventh	AGHS2	AGHS3	1.3087	0.4861	73.62	2.69	0.0088	0.0893
Year Seventh	AGHS2	AGHS4	0.9115	0.4861	73.62	1.88	0.0647	0.4256
Year Seventh	AGHS3	AGHS4	-0.3972	0.4861	73.62	-0.82	0.4165	0.9635

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Family	Estimate
Year Seventh	AGHS2	5.2601	A
Year Seventh			A
Year Seventh	AGHS1	4.8530	A
Year Seventh			A
Year Seventh	AGHS4	4.3486	A
Year Seventh			A
Year Seventh	AGHS3	3.9514	A
Year Seventh
Year Seventh	AGH2	2.2175	B
Year Seventh			B
Year Seventh	AGH25	2.0697	B

Year*Family F Test

F Test for Year*Family Least Squares Means Slice
Slice	Num DF	Den DF	F Value	Pr > F
Year Sixth	5	73.62	13.84	<.0001

Year*Family Diffs

Simple Differences of Year*Family Least Squares Means Adjustment for Multiple Comparisons: Tukey
Slice	Family	_Family	Estimate	Standard Error	DF	t Value	Pr > |t|	Adj P
Year Sixth	AGH2	AGH25	-0.4612	0.4861	73.62	-0.95	0.3459	0.9322
Year Sixth	AGH2	AGHS1	-2.7922	0.4861	73.62	-5.74	<.0001	<.0001
Year Sixth	AGH2	AGHS2	-3.1903	0.4861	73.62	-6.56	<.0001	<.0001
Year Sixth	AGH2	AGHS3	-2.1436	0.4861	73.62	-4.41	<.0001	0.0005
Year Sixth	AGH2	AGHS4	-2.1314	0.4861	73.62	-4.38	<.0001	0.0005
Year Sixth	AGH25	AGHS1	-2.3311	0.4861	73.62	-4.80	<.0001	0.0001
Year Sixth	AGH25	AGHS2	-2.7291	0.4861	73.62	-5.61	<.0001	<.0001
Year Sixth	AGH25	AGHS3	-1.6824	0.4861	73.62	-3.46	0.0009	0.0113
Year Sixth	AGH25	AGHS4	-1.6703	0.4861	73.62	-3.44	0.0010	0.0122
Year Sixth	AGHS1	AGHS2	-0.3981	0.4861	73.62	-0.82	0.4155	0.9631
Year Sixth	AGHS1	AGHS3	0.6486	0.4861	73.62	1.33	0.1862	0.7650
Year Sixth	AGHS1	AGHS4	0.6608	0.4861	73.62	1.36	0.1782	0.7507
Year Sixth	AGHS2	AGHS3	1.0467	0.4861	73.62	2.15	0.0346	0.2725
Year Sixth	AGHS2	AGHS4	1.0589	0.4861	73.62	2.18	0.0326	0.2606
Year Sixth	AGHS3	AGHS4	0.01218	0.4861	73.62	0.03	0.9801	1.0000

Year*Family Lines

Tukey Grouping for Year*Family Least Squares Means Slice (Alpha=0.05)
LS-means with the same letter are not significantly different.
Slice	Family	Estimate
Year Sixth	AGHS2	4.4584	A
Year Sixth			A
Year Sixth	AGHS1	4.0603	A
Year Sixth			A
Year Sixth	AGHS3	3.4117	A
Year Sixth			A
Year Sixth	AGHS4	3.3995	A
Year Sixth
Year Sixth	AGH25	1.7292	B
Year Sixth			B
Year Sixth	AGH2	1.2681	B

Survival rates

Sur1 <- ggplot(hane_survival_year,
               aes(x = Survival, y = Rate, fill = Survival, Color = Survival)) +
  stat_summary(fun = mean, geom = "bar", position = "dodge") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 100), breaks = c(0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("green4", "red3"), 0.5)) +
  labs(title = "Survival Rate by Growing Season") +
  xlab("") +
  ylab("Survival rate (%)") +
  facet_wrap( ~ Year) +
  guides(fill = "none") +
  theme_classic()

plotly::ggplotly(Sur1)

Sur2 <- ggplot(hane_survival_year_family,
               aes(x = Family, y = Rate, fill = Survival, color = Survival)) +
  stat_summary(fun = mean, geom = "bar", position = "dodge") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 100), breaks = c(0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("green4", "red3"), 0.5)) +
  labs(title = "Survival Rate by Growing Season and Family") +
  xlab("") +
  ylab("Survival rate (%)") +
  facet_wrap( ~ Year) +
  guides(fill = "none") +
  theme_classic()

plotly::ggplotly(Sur2)

Sur3 <- ggplot(hane_survival_year_location,
               aes(x = Location, y = Rate, fill = Survival, Color = Survival)) +
  stat_summary(fun = mean, geom = "bar", position = "dodge") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 100), breaks = c(0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("green4", "red3"), 0.5)) +
  labs(title = "Survival Rate by Growing Season and Location") +
  xlab("") +
  ylab("Survival rate (%)") +
  facet_wrap( ~ Year) +
  guides(fill = "none") +
  theme_classic()

plotly::ggplotly(Sur3)

Sur4 <- ggplot(hane_survival_year_location_family,
               aes(x = Family, y = Rate, fill = Survival, color = Survival)) +
  stat_summary(fun = mean, geom = "bar", position = "dodge") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 100), breaks = c(0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("green4", "red3"), 0.5)) +
  labs(title = "Survival Rate by Growing Season, Location, and Family") +
  xlab("") +
  ylab("Survival rate (%)") +
  facet_wrap( ~ ID, ncol = 3) +
  guides(fill = "none") +
  theme_classic()

plotly::ggplotly(Sur4)

GLD

Fig1 <- ggplot(hane_gld,
                 aes(x = Family, y = GLD, fill = Family, color = Family)) +
  stat_summary(fun.data = mean_sd, geom = "errorbar", width = 0.2, color = "grey") +
  stat_summary(fun = mean, geom = "bar") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 60), breaks = c(0, 10, 20, 30, 40, 50, 60), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("#E69F00", "#56B4E9", "#009E73",
                                     "#F0E442", "#D55E00", "#CC79A7"), 0.5)) +
  labs(title = "GLD: Family Effect") +
  xlab("Family") +
  ylab("GLD (mm)") +
  annotate("text",
           x = c(1, 2, 3, 4, 5, 6),
           y = c(50),
           label = c("b", "ab", "a", "a", "ab", "ab")) +
  guides(fill = "none") +
  theme_classic()

plotly::ggplotly(Fig1)

Fig2 <- ggplot(hane_gld,
               aes(x = Year, y = GLD, fill = Location, color = Location)) +
  stat_summary(fun.data = mean_sd, geom = "errorbar", position = position_dodge(0.9), width = 0.2, color = "grey") +
  stat_summary(fun = mean, geom = "bar", position = "dodge") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 70), breaks = c(0, 10, 20, 30, 40, 50, 60, 70), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("#000000", "#999999"), 0.5)) +
  labs(title = "GLD: Year * Location") +
  xlab("Growing Season") +
  ylab("GLD (mm)") +
  annotate("text",
           x = c(1.78, 2.22, 2.78, 3.22),
           y = c(20, 20, 55, 50),
           label = c("b", "a", "a", "b")) +
  guides(fill = "none") +
  theme_classic()

plotly::ggplotly(Fig2)

Fig3 <- ggplot(hane_gld,
               aes(x = Year, y = GLD, fill = Family, color = Family)) +
  stat_summary(fun.data = mean_sd, geom = "errorbar", position = position_dodge(0.9), width = 0.2, color = "grey") +
  stat_summary(fun = mean, geom = "bar", position = "dodge") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 80), breaks = c(0, 10, 20, 30, 40, 50, 60, 70, 80), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("#E69F00", "#56B4E9", "#009E73",
                               "#F0E442", "#D55E00", "#CC79A7"), 0.5)) +
  labs(title = "GLD: Year * Family") +
  xlab("Growing Season") +
  ylab("GLD (mm)") +
  annotate("text",
           x = c(2.62, 2.78, 2.93, 3.08, 3.22, 3.38),
           y = c(45, 50, 55, 60, 55, 55),
           label = c("b", "b", "a", "a", "a", "a")) +
  guides(fill = "none") +
  theme_classic()

plotly::ggplotly(Fig3)

Height

Fig4 <- ggplot(hane_ht,
               aes(x = Family, y = Height, fill = Family, color = Family)) +
  stat_summary(fun.data = mean_sd, geom = "errorbar", width = 0.2, color = "grey") +
  stat_summary(fun = mean, geom = "bar") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 12), breaks = c(0, 2, 4, 6, 8, 10, 12), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("#E69F00", "#56B4E9", "#009E73",
                                     "#F0E442", "#D55E00", "#CC79A7"), 0.5)) +
  labs(title = "Height: Year * Location") +
  xlab("Family") +
  ylab("Height (m)") +
  annotate("text",
           x = c(1, 2, 3, 4, 5, 6),
           y = c(10),
           label = c("d", "cd", "a", "a", "bc", "ab")) +
  guides(fill = "none") +
  theme_classic()

plotly::ggplotly(Fig4)

Fig5 <- ggplot(hane_ht,
                  aes(x = Year, y = Height, fill = Location, color = Location)) +
  stat_summary(fun.data = mean_sd, geom = "errorbar", position = position_dodge(0.9), width = 0.2, color = "grey") +
  stat_summary(fun = mean, geom = "bar", position = "dodge") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 12), breaks = c(0, 2, 4, 6, 8, 10, 12), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("#000000", "#999999"), 0.5)) +
  labs(title = "Height: Year * Location" ) +
  xlab("Growing Season") +
  ylab("Height (m)") +
  annotate("text",
           x = c(3.78, 4.22, 4.78, 5.22, 5.78, 6.22, 6.78, 7.22),
           y = c(5, 4, 7, 5.5, 8.5, 7.5, 11, 7.5),
           label = c("a", "b", "a", "b", "a", "b", "a", "b")) +
  theme_classic()

plotly::ggplotly(Fig5)

Fig6 <- ggplot(hane_ht,
               aes(x = Year, y = Height, fill = Family, color = Family)) +
  stat_summary(fun.data = mean_sd, geom = "errorbar", position = position_dodge(0.9), width = 0.2, color = "grey") +
  stat_summary(fun = mean, geom = "bar", position = "dodge") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 12), breaks = c(0, 2, 4, 6, 8, 10, 12), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("#E69F00", "#56B4E9", "#009E73",
                                     "#F0E442", "#D55E00", "#CC79A7"), 0.5)) +
  labs(title = "Height: Year * Family") +
  xlab("Growing Season") +
  ylab("Height (m)") +
  annotate("text",
           x = c(2.63, 2.77, 2.92, 3.06, 3.22, 3.38,
                 3.63, 3.77, 3.92, 4.06, 4.22, 4.38,
                 4.63, 4.77, 4.92, 5.06, 5.22, 5.38,
                 5.63, 5.77, 5.92, 6.06, 6.22, 6.38,
                 6.63, 6.77, 6.92, 7.06, 7.22, 7.38,
                 7.63, 7.77, 7.92, 8.06, 8.22, 8.38
                 ),
           y = c(2.2, 2.5, 3, 3.3, 2.7, 2.9,
                 3.1, 3.3, 4.3, 4.5, 3.7, 3.9,
                 5, 5.5, 6.5, 6.8, 5.9, 6.2,
                 6.5, 7, 8.1, 8.4, 7.4, 7.8,
                 7.7, 8.5, 10.4, 10.4, 9.8, 10.2,
                 9.3, 9.2, 10.5, 10.4, 10, 10.3
                 ),
           label = c("c", "bc", "a", "a", "abc", "ab",
                     "c", "bc", "a", "a", "bc", "ab",
                     "d", "cd", "ab", "a", "bc", "ab",
                     "c", "bc", "a", "a", "b", "a",
                     "c", "c", "a", "a", "b", "ab",
                     "dc", "d", "a", "a", "bc", "ab"
                     )) +
  theme_classic()

plotly::ggplotly(Fig6)

DBH

Fig7 <- ggplot(hane_dbh,
               aes(x = Family, y = DBH, fill = Family, color = Family)) +
  stat_summary(fun.data = mean_sd, geom = "errorbar", width = 0.2, color = "grey") +
  stat_summary(fun = mean, geom = "bar") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 20), breaks = c(0, 5, 10, 15, 20), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("#E69F00", "#56B4E9", "#009E73",
                                     "#F0E442", "#D55E00", "#CC79A7"), 0.5)) +
  labs(title = "DBH: Family Effect") +
  xlab("Family") +
  ylab("DBH (cm)") +
  annotate("text",
           x = c(1, 2, 3, 4, 5, 6),
           y = c(15),
           label = c("d", "cd", "b", "b", "cb", "a")) +
  theme_classic()

plotly::ggplotly(Fig7)

Fig8 <- ggplot(hane_dbh,
               aes(x = Year, y = DBH, fill = Location, color = Location)) +
  stat_summary(fun.data = mean_sd, geom = "errorbar", position = position_dodge(0.9), width = 0.2, color = "grey") +
  stat_summary(fun = mean, geom = "bar", position = "dodge") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 20), breaks = c(0, 5, 10, 15, 20), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("#000000", "#999999"), 0.5)) +
  labs(title = "DBH: Year * Location" ) +
  xlab("Growing Season") +
  ylab("DBH (cm)") +
  annotate("text",
           x = c(0.78, 1.22, 1.78, 2.22, 2.78, 3.22, 3.78, 4.22, 4.78, 5.22),
           y = c(5, 4, 8, 7, 11, 10, 13, 12, 14, 13),
           label = c("a", "b", "a", "b", "a", "b", "a", "b", "a", "b")) +
  guides(fill = "none") +
  theme_classic()

plotly::ggplotly(Fig8)

Fig9 <- ggplot(hane_dbh,
               aes(x = Year, y = DBH, fill = Family, color = Family)) +
  stat_summary(fun.data = mean_sd, geom = "errorbar", position = position_dodge(0.9), width = 0.2, color = "grey") +
  stat_summary(fun = mean, geom = "bar", position = "dodge") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 20), breaks = c(0, 5, 10, 15, 20), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("#E69F00", "#56B4E9", "#009E73",
                                     "#F0E442", "#D55E00", "#CC79A7"), 0.5)) +
  labs(title = "DBH: Year * Family") +
  xlab("Growing Season") +
  ylab("DBH (cm)") +
  annotate("text",
           x = c(0.63, 0.77, 0.92, 1.06, 1.22, 1.38,
                 1.63, 1.77, 1.92, 2.06, 2.22, 2.38,
                 2.63, 2.77, 2.92, 3.06, 3.22, 3.38,
                 3.63, 3.77, 3.92, 4.06, 4.22, 4.38,
                 4.63, 4.77, 4.92, 5.06, 5.22, 5.38
           ),
           y = c(3.2, 3.5, 5, 5.3, 4.7, 4.9,
                 6.1, 6.3, 7.8, 7.9, 7.5, 8.5,
                 8.9, 9.3, 10.1, 10.2, 10.8, 11.8,
                 11.1, 11.3, 11.9, 12.3, 12.8, 13.2,
                 12.7, 12.7, 12.9, 13.1, 13.5, 14.7
           ),
           label = c("c", "c", "ab", "a", "bc", "a",
                     "d", "cd", "b", "b", "bc", "a",
                     "c", "c", "b", "b", "b", "a",
                     "d", "cd", "bcd", "bc", "b", "a",
                     "b", "b", "b", "b", "b", "a"
           )) +
  guides(fill = "none") +
  theme_classic()

plotly::ggplotly(Fig9)

Crown

Fig10 <- ggplot(hane_crown,
                aes(x = Family, y = Crown, fill = Family, color = Family)) +
  stat_summary(fun.data = mean_sd, geom = "errorbar", width = 0.2, color = "grey") +
  stat_summary(fun = mean, geom = "bar") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 10), breaks = c(0, 2, 4, 6, 8, 10), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("#E69F00", "#56B4E9", "#009E73",
                                     "#F0E442", "#D55E00", "#CC79A7"), 0.5)) +
  labs(title = "Crown Height: Family Effect") +
  xlab("Family") +
  ylab("Crown Height (m)") +
  annotate("text",
           x = c(1, 2, 3, 4, 5, 6),
           y = c(8),
           label = c("b", "b", "a", "a", "a", "a")) +
  guides(fill = "none") +
  theme_classic()

plotly::ggplotly(Fig10)

Fig11 <- ggplot(hane_crown,
               aes(x = Year, y = Crown, fill = Location, color = Location)) +
  stat_summary(fun.data = mean_sd, geom = "errorbar", position = position_dodge(0.9), width = 0.2, color = "grey") +
  stat_summary(fun = mean, geom = "bar", position = "dodge") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 10), breaks = c(0, 2, 4, 6, 8, 10), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("#000000", "#999999"), 0.5)) +
  labs(title = "Crown Height: Year * Location" ) +
  xlab("Growing Season") +
  ylab("Crown Height (m)") +
  annotate("text",
           x = c(0.78, 1.22, 1.78, 2.22, 2.78, 3.22),
           y = c(3.2, 1.7, 6.5, 4, 7.3, 5),
           label = c("a", "b", "a", "b", "a", "b")) +
  guides(fill = "none") +
  theme_classic()

plotly::ggplotly(Fig11)

Fig12 <- ggplot(hane_crown,
               aes(x = Year, y = Crown, fill = Family, color = Family)) +
  stat_summary(fun.data = mean_sd, geom = "errorbar", position = position_dodge(0.9), width = 0.2, color = "grey") +
  stat_summary(fun = mean, geom = "bar", position = "dodge") +
  font("title", size = 22) +
  font("xlab", size = 18) +
  font("ylab", size = 18) +
  font("xy.text", size = 16) +
  scale_y_continuous(limits = c(0, 10), breaks = c(0, 2, 4, 6, 8, 10), expand = c(0,0)) +
  scale_fill_manual(values = alpha(c("#E69F00", "#56B4E9", "#009E73",
                                     "#F0E442", "#D55E00", "#CC79A7"), 0.5)) +
  labs(title = "Crown Height: Year * Family") +
  xlab("Growing Season") +
  ylab("Crown Height (m)") +
  annotate("text",
           x = c(1.63, 1.77, 1.92, 2.06, 2.22, 2.36,
                 2.63, 2.77, 2.92, 3.06, 3.22, 3.36
           ),
           y = c(2.5, 3.4, 6, 6.5, 5.5, 5.7,
                 4, 3.9, 7.2, 7.3, 6.1, 6.6
           ),
           label = c("b", "b", "a", "a", "a", "a",
                     "b", "b", "a", "a", "a", "a"
           )) +
  guides(fill = "none") +
  theme_classic()

plotly::ggplotly(Fig12)